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Section  A-l 
INTRODUCTION 

A  substantial  portion  of  the  Ship  Identification  Study  (SIS)  effort 
was  devoted  to  an  investigation  of  atmospheric  optical  effects  likely  to  be 
encountered  in  the  long  slant  paths  required  for  ship  identification  at 
acceptable  stand-off  ranges.  Potentially  troublesome  atmospheric  effects 
can  occur  to  a  laser  beam  utilized  for  illumination  in  an  active  system,  and 
to  reflected  energy  returned  to  the  observation  aircraft  from  the  target  area. 
(In  some  respects  the  significant  effects  for  the  illuminating  beam  are  dif¬ 
ferent  from  those  for  the  returned  energy.) 

There  were  three  principal  reasons  for  the  large  effort  expended  in  the 
investigation  of  atmospheric  phenomena  in  the  optical  wavelength*  region: 

(1)  It  was  anticipated  that  atmospheric  phenomena  would  play  a  large 
role  in  determining  the  potential  merits  of  the  many  candidate  sensor  concepts 
to  be  considered; 

(2)  Atmospheric  molecular  absorption  coefficients  for  laser  frequencies 
in  a  potentially  attractive  portion  of  the  optical  spectrum  (3.5  to  4.0  micro¬ 
meter)  were  not  well-established,  to  the  extent  that  any  calculated  estimate 
was  considered  suspect  and  controversial;  and 

(3)  A  potentially  troublesome  effect,  optical  turbulence  (both  ambient 
and  aircraf t- induced) ,  was  but  poorly  known  and  little  understood  for  long 
air-to-ground  slant  paths. 

It  was  known  that  the  problems  indicated  by  (2)  and  (3)  above  were 
subjects  of  rather  intensive  research  at  several  laboratories.  Researchers 
at  those  laboratories  were  consulted  during  the  course  of  the  study,  and  very 
significant  advance  unpublished  results  were  thereby  obtained  for  application 
to  the  SIS  effort.  Thus  the  findings  of  this  report  depended  heavily  upon 


*As  used  herein,  the  term  "optical  wavelengths"  encompasses  electromagnetic 
waves  having  free-space  wavelengths  in  the  region  from  0.3  to  15.0  micrometers. 
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information  supplied  by  the  Optical  Physics  Laboratory  of  the  Air  Force  Cam¬ 
bridge  Research  Laboratories,  the  Optical  Sciences  Division  of  the  Naval 
Research  Laboratory,  Air  Force  Weapons  Laboratory  divisions  under  the  Advanced 
Radiation  Technology  Office,  the  Electro  Science  Laboratory  of  Ohio  State 
University,  and  Science  Applications  Incorporated.  (Information  from  the 
latter  two  organizations  resulted  from  efforts  sponsored  by  Defense  Research 
Projects  Agency  under  Rome  Air  Development  Center  contracts,  and  release  of 
advance  data  from  those  efforts  was  arranged  through  RADC.)  Our  thanks  to  the 
many  individuals  involved  who  shared  their  expertise  with  us  for  this  project. 


References  1  through  10  of  this  appendix  are  formal  published  reports 
in  the  open  literature  which  constitute  an  excellent  starting  point  and  pack¬ 
age  of  data  sources  for  beginning  a  study  of  atmospheric  effects  upon  an 
electro-optical  (EO)  system.  These  references  were  so  employed  in  the  SIS 
effort.  A  few  papers  from  technical  journals  also  proved  helpful,  and  these 
are  listed.  A  characteristic  of  this  study,  however,  is  the  large  use  that 
was  made  of  advance  unpublished  data  obtained  from  researchers  in  the  field. 
Many  of  the  listed  references  relate  to  such  data.  A  bibliography  is  pro¬ 
vided  for  sources  not  cited  specifically. 


As  various  aspects  of  the  atmospheric  transmission  problem  were  investi¬ 
gated,  draft  papers  were  written  to  document  the  individual  efforts.  Sections 
A-3  through  A-ll  of  this  appendix  consist  essentially  of  those  draft  papers, 
with  some  revision  and  updating  where  appropriate. 


A  "66-kilometer  slant  path"  is  referred  to  throughout  this  appendix 
in  the  context  of  a  nominal  atmospheric  path  over  which  SIS  must  be  applicable. 
This  nominal  slant  path  is  defined  in  the  sketch  below. 
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Section  A-2 

ATMOSPHERIC  EFFECTS  SYNOPSIS 


GENERAL 

Optical  energy  reflected  or  emitted  from  a  target  element  and  received 
by  a  sensor  element  at  a  distant  location  is  affected  in  numerous  ways  by 
transmission  through  the  atmosphere.  Phenomena  of  major  importance  in  the 
visible  and/or  IR  wavelength  regions  are  (1)  attenuation  of  the  target  energy 
by  path  absorption  and  scattering  out  of  the  path-,  (2)  scattering  of  ex¬ 
traneous  background  energy  into  the  transmission  path,  (3)  thermal  radiance  of 
the  air  mass  in  the  path,  and  (4)  optical  turbulence  within  the  path  air  mass. 
For  a  collimated  or  focused  laser  beam  used  for  target  illumination,  (1)  and 
(4)  are  the  atmospheric  phenomena  of  major  importance. 

Atmospheric  absorption,  scattering,  and  thermal  radiance  may  be  cal¬ 
culated  approximately  by  the  use  of  appropriate  atmospheric  models,  together 
with  mathematical  models  representing  the  atmospheric  effects  phenomena.  SRL 
mainly  uses  models  and  coefficient  tables  developed  by  the  Optical  Physics 
Laboratory  of  AFCRL,  as  documented  in  several  of  the  cited  references.  We 
also  employ  computer  programs  and  data  tables  developed  by  AFCRL  for  computer 
calculation  of  spectral  band  atmospheric  transmittance,  manual  calculation  of 
which  is  tedious  and  costly. 

The  several  scattering  and  absorption  coefficients  are  strong  functions 
of  altitude,  and  this  factor  must  be  accounted  for  in  an  appropriate  fashion 
in  transmittance,  path  radiance,  and  optical  turbulence  calculations  for  air- 
to-surface  slant  paths. 

The  effects  of  some  adverse  weather  conditions  (e.g.,  rain  and  fog)  are 
predictable  to  a  limited  extent  if  the  relevant  parameters  (e.g.,  rainfall 
rate  and  extent)  can  be  measured  or  estimated.  A  USAF  Project  RAND  report, 
R-1523-PR  (Ref.  11),  provides  information  for  calculating  attenuation  through 
rain,  and  the  calculated  results  appear  to  be  in  reasonable  agreement  with 


measured  values.  It  should  be  noted  that  rain  attenuation  calculations  based 
upon  the  approach  in  Middleton's  classic  text  (Ref.  12)  give  totally  erro¬ 
neous  results,  and  should  not  be  used. 

The  McClatchey  models  (AFCRL)  (Refs  1-8)  do  not  take  path  optical  tur¬ 
bulence  into  account,  and  this  phenomenon  must  be  accounted  for  separately  to 
obtain  quantitative  estimates  of  the  effects  of  this  factor.  The  theoretical 
basis  in  this  area  is  in  a  primitive  state  of  development;  however,  quasi- 
empirical  mathematical  models  are  available  for  making  engineering  estimates, 
and  such  computed  estimates  appear  to  agree  reasonably  well  with  experimental 
results  over  a  nominal  range  of  atmospheric  conditions  and  target-to- 
receiver  geometries. 

Extensive  use  was  made  of  a  recent  RAND  Corporation  work  (Ref.  13)  in 
this  area,  supplemented  by  basic  data  from  Tatarski  (Ref.  14)  and  other 
researchers.  Experimental  data  relating  to  ambient  altitude  variations  and 
aircraft- induced  turbulence  were  obtained  from  the  AFWL  Airborne  Laser 
Laboratory  program  and  applied  to  the  SIS  investigation.  More  data  are 
expected  to  be  available  from  AFWL  during  1976. 


VISIBLE  AND  NEAR- I R  WAVELENGTH  REGION 

In  this  wavelength  region,  optical  attenuation  and  path  radiance  are 
mainly  attributable  to  scattering  of  energy  by  air  molecules  and  by  atmos¬ 
pheric  haze  caused  by  suspended  dust,  smoke  and  other  foreign  matter  (collec¬ 
tively  termed  "aerosols"). 

At  the  shorter  wavelengths,  Rayleigh  scattering  (by  air  molecules)  and 
Mie  scattering  (by  atmospheric  aerosols)  are  of  somewhat  comparable  signifi¬ 
cance  in  the  lower  atmosphere  on  a  relatively  clear  day.  On  a  relatively  hazy 
day,  aerosol  scattering  greatly  predominates  and  becomes  the  limiting  factor 
in  visible  wavelength  transmission.  Aerosol  absorption  is  marginally 
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significant  (but  is  often  ignored),  while  molecular  absorption  is  negligible, 
except  as  noted  in  the  following  paragraph,  and  may  be  ignored  for  band  cal¬ 
culations  in  the  visible  wavelength  region. 

Water  vapor  molecules  are  not  important  scatterers  or  absorbers  of 
visible  wavelengths;  hence,  visible  transmittance  is  not  affected  signifi¬ 
cantly  by  atmospheric  humidity.  That  situation  changes  at  the  red  end  of 
the  visible  band,  however,  with  significant  water  vapor  absorption  starting 
at  about  0.69  micrometer,  as  illustrated  in  Sets  I  and  IV  of  the  trans¬ 
mittance  graphs  in  Section  A-6.  Those  graphs  show  computer-calculated  band 
transmittance  over  an  SIS  air-to-ground  slant  path,  for  the  visible  to  near- 
IR  wavelength  region  of  0.3  to  1.8  micrometers,  for  a  wide  range  of  atmospheric 
conditions . 

With  daylight  illumination,  path  radiance  at  visible  wavelengths 
consists  mainly  of  sky  background  energy  which  is  scattered  into  the  sensor 
field  of  view  (FOVl.  With  artificial  illumination  in  darkness,  backscattered 
energy  from  the  illuminator  is  the  principal  source  of  path  radiance. 

In  this  wavelength  region  optical  turbulence  can  cause  a  substantial 
increase  in  the  imaging  point  spread  function,  with  a  concomitant  degradation 
in  the  image  quality  achievable  in  high-resolution  photography  and  television. 
For  active  night-time  systems  employing  laser  illuminators,  optical  turbulence 
induces  beam  wander  and  beam  spread  into  the  collimated  laser  beam,  effectively 
placing  a  minimum  diameter  limit  upon  the  useful  size  (and  associated  in¬ 
tensity)  of  the  beam.  Optical  turbulence  effects  vary  over  a  wide  range  in 
the  lower  atmosphere,  depending  upon  the  thermal  dynamics  of  the  air  and 
adjacent  surfaces. 

8-13  MICROMETER  IR  BAND 

In  this  wavelength  region  optical  attenuation  is  chiefly  attributable 
to  molecular  absorption  by  carbon  dioxide  and  water  vapor,  and  path  radiance 
is  mainly  due  to  thermal  emission  of  the  air  mass  in  the  transmission  path. 
Absorption  and  scattering  by  atmospheric  aerosols  are  of  some  significance 


but  are  minor  contributors  to  the  total  except  under  conditions  of  rel  itively 
heavy  haze  and/or  low  atmospheric  temperature  and  water  vapor  content. 


Under  most  atmospheric  conditions  the  predominant  factor  influencing 
transmittance  in  this  wavelength  region  is  atmospheric  water  vapor  content 
(absolute,  not  relative,  humidity).  Also,  since  water  vapor  absorption  is 
strongly  temperature  dependent,  the  air  temperature  of  the  transmission  path 
is  a  significant  secondary  factor.  (A  warm,  humid  atmosphere  provides 
poorest  transmittance.) 

In  this  wavelength  region  optical  turbulence  appears  to  be  much  less 
significant  than  at  visible  wavelengths.  Minor  degradations  in  the  imaging 
point  spread  function  and  illumination  beam  parameters  may  occur  under 
conditions  of  moderate  to  strong  turbulence. 

3-5  MICROMETER  IR  BAND 

In  this  wavelength  region  the  atmosphere  exhibits  optical  properties 
which  are  a  composite  of  those  previously  cited  for  shorter  and  longer  wave¬ 
length  bands.  Molecular  and  aerosol  absorption  and  aerosol  scattering  are  all 
significant  in  the  attenuation  equation,  and  both  scattering  and  air  mass 
thermal  radiation  should  be  considered  in  estimating  path  radiance.  The 
effects  of  atmospheric  optical  turbulence  upon  imaging  and  laser  beam  trans¬ 
mission  are  significantly  less  than  they  are  in  the  visible  region,  but 
somewhat  greater  than  they  are  in  the  8-13  micrometer  band. 

Portions  of  this  band  (in  the  vicinity  of  2.2  and  3.8  micrometers) 
exhibit  optimum  optical  transmission  over  a  wide  range  of  atmospheric  variables. 
Optical  energy  in  these  wavelength  regions  penetrates  atmospheric  haze  much 
better  than  shorter  (e.g.,  visible)  wavelengths,  and  is  much  less  affected  by 
atmospheric  water  vapor  than  longer  wavelengths  (e.g.,  the  8-13  micrometer 
IR  region) . 
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ATMOSPHERIC  IMPACT  UPON  LASER  BEAM  ILLUMINATORS 

Atmospheric  attenuation  directly  decreases  the  illumination  energy 
available  in  the  target  region.  Atmospheric  turbulence,  whether  aircraft- 
induced  or  ambient,  breaks  up  the  smoothly-varying  spatial  pattern  of  a 
collimated  or  focused  beam  and  produces  "scintillation"  or  rapidly  varying 
intensity  fluctuations  within  the  beam.  A  random  wandering  of  the  beam  about 
a  central  point,  together  with  skew,  astigmatism,  and  spreading,  are 
additional  beam  distortions  produced  by  optical  turbulence.  The  end  result 
is  a  spread,  wandering,  fluctuating  beam  whose  average  intensity  on  target 
can  be  much  reduced  compared  to  that  predicted  by  diffraction  limit  theory 
for  vacuum  transmissions.  This  degradation  in  beam  quality  can  be  very 
significant  over  long  atmospheric  paths  (particularly  if  sharp,  efficient 
illuminating  beams  are  desired)  and  must  be  accounted  for  in  system  design. 

ATMOSPHERIC  IMPACT  UPON  SENSORS 

Atmospheric  attenuation  decreases  the  optical  signal  at  the  receiver, 
while  path  radiance  introduces  a  background  "noise"  level.  These  factors 
directly  influence  the  available  contrast  in  a  sensor  (such  as  photographic 
film)  which  detects  total  energy  content.  In  an  electro-optical  sensor 
(such  as  TV  or  FLIR)  the  path  radiance  can  be  largely  suppressed  by  filtering, 
with  mainly  the  video  or  ac  component  recorded  or  displayed.  Scene  contrast 
at  the  sensor  (and  signal-to-noise  ratio)  are  then  influenced  mainly  by  atmos¬ 
pheric  attenuation  and  sensor  detectivity.  Atmospheric  turbulence  over  the 
transmission  path  can  be  expected  to  produce  radiometer/photometer  reading 
fluctuations  and  noticeable  TV  or  FLIR  image  distortion  and  resolution  degrada¬ 
tion  some  significant  percentage  of  the  time. 

SOME  OBSERVATIONS  ABOUT  CALCULATION  ACCURACIES  AND  SCALING 

-A 

The  molecular  scattering  coefficient,  a  ,  exhibits  a  X  dependence 
upon  wavelength,  and  the  coefficient  scales  linearly  with  the  number  of 
molecules  in  the  path  (and  hence  approximately  with  air  density).  Thus 
attenuation  attributable  to  molecular  scattering  is  easily  and  accurately 
calculable. 


For  a  given  constituency  of  materials  and  particle  sizes  making  up 
an  atmospheric  aerosol  model,  the  aerosol  scattering  coefficient,  a  ,  and 
aerosol  absorption  coefficient,  k  ,  scale  linearly  with  the  number  of 

cl 

particles  in  the  path,  and  the  variation  with  wavelength  is  reasonably  well 
behaved.  Hence,  aerosol  scattering  and  attenuation  are  easily  calculated 
for  model  atmospheres  for  which  a  and  k  tables  have  been  established  for 

cl  3 

a  specified  constituency  and  concentration  of  aerosol  particles  (including 
variation  with  altitude).  The  accuracy  of  particular  atmospheric  aerosol 
path  transmittance  estimates  based  upon  such  calculations  is  likely  to  be 
rather  poor,  however,  because  the  actual  aerosol  constituency  and  dis¬ 
tributions  present  are  not  likely  to  be  known  very  well  and  hence  will  not 
fit  the  model  very  well.  This  factor  is  probably  the  principal  source  of 
error  in  estimating  long  path  atmospheric  effects  at  visible  wavelengths. 

The  molecular  absorption  coefficient,  km,  is  a  complicated  function 
of  the  total  pressure,  partial  pressures,  temperature,  and  specific  con¬ 
stituency  of  the  atmospheric  gases  in  the  transmission  path,  and  of  optical 
wavelength.  Calculation  of  this  factor  in  atmospheric  transmittance  is 
difficult  and  often  the  results  are,  at  best,  rough  estimates.  A  published 
atmospheric  absorption  coefficient  for  one  temperature,  pressure,  and  con¬ 
stituency  is  not  simply  scalable  to  another  set  of  conditions  because 
individual  air  molecule  constituents  scale  differently  with  partial  pressure, 
total  pressure,  and  temperature,  while  wavelength  dependence  is  sharply 
structured  and  relatively  unpredictable  (except  by  machine  computation). 
Figures  A-l  through  A-4  (from  Ref.  15)  illugtrate  the  complexity  of  this 
factor  in  atmospheric  transmittance.  See  also,  typical  high-resolution  spec¬ 
tral  transmittance  curves  in  the  3.8  micrometer  region  (Figure  A- 22  of  Section 
A-7)  and  in  the  10.6  micrometer  region  (Figure  A- 31  of  Section  A-8) . 

It  is  because  of  the  problem  indicated  in  the  foregoing  paragraph,  that 
a  large  body  of  spectral  line  absorptior  coefficients  and  a  band  transmittance 
model  (LOWTRAN)  for  a  family  of  model  atmospheres  have  been  developed  by 
AFCRL.  Engineering  estimates  for  practical  application  purposes  can  be 
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Figure  A-l.  Contributor,  to  the  Molecular  Absorption  of  the  P  (9) 
DF  Line  (Mid-Latitude  Summer,  Sea  Level)  1 
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Figure  A-2.  Temperature  Dependence  of  the  Major  Molecular  Line 
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computed  for  "typical"  atmospheric  conditions  over  a  rather  wide  range, 
by  application  of  the  AFCRL  models  and  data  (using  either  manual  or  machine 
computation) . 


DOPPLER  SHIFT 

The  position  of  certain  laser  lines  r,  lative  to  high-resolution 
atmospheric  transmittance  spectra  (e.g.,  Figure  A-4)  gave  rise  to  this  ques¬ 
tion:  Is  the  doppler  shift  in  frequency  significant  with  respect  to  atmos¬ 
pheric  attenuation  of  a  laser  beam  transmitted  from  an  airplane  in  flight? 

A  simple  calculation  showed  that  the  maximum  doppler  shift  (transmit 
direction  at  0°  or  180°  to  aircraft  velocity)  for  an  aircraft  velocity  of 
500  knots  is: 

Af  =  8.574  x  10" 7  f 

o 

where  f  is  laser  frequency.  For  the  DF  2-1  P7  line  (for  example)  the 
doppler  shift  at  500  knots  is: 

Af  =  8.574  x  10”7  x  2655.863  per  cm 

=  0.00228  waves/cm 

This  is  less  than  the  order  of  accuracy  (0.003/cm)  stated  for  the  best 
available  DF  laser  frequency  measurements  (Ref.  16).  Comparing  this  value 
with  transmittance  spectra  available  from  OSU,  AFCRL,  and  SAI,  it  is  con¬ 
cluded  that  the  effect  would  be  a  marginal  one,  of  no  practical  significance 
with  respect  to  transmittance  from  a  B-52.  It  appears  that  the  shift  would 
start  to  be  significant  (for  a  few  DF  laser  lines)  at  about  0.01  cm  ^  Af,  or 
a  velocity  of  about  2000  knots.  Doppler  shift  would  probably  cause  a  small 
(but  insignificant)  improvement  in  the  transmittance  of  most  CO2  laser  lines, 
because  the  shift  would  be  away  from  the  peak  of  the  related  CO2  absorption 
line  in  the  atmosphere. 
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Section  A- 3 

PATH  RADIANCE  AND  TRANSMITTANCE 

PATH  EQUATION  FOR  RECEIVED  POWER 

Assume  that  emitted  or  reflected  energy  from  a  target  element  of  uni¬ 
form  radiance  is  focused'  upon  a  sensor  surface  or  detector  element  of  a  photom¬ 
eter,  radiometer,  or  other  receiver  which  is  at  some  distance  from  the  target. 
For  a  horizontal  path  of  uniform  characteristics,  the  power  received  by  the 
sensor  element  can  be  written  as  follows: 


PR  =  N  T  (K)  +  N,  (SF)  (K)  +  N  Em  (K) 
K  t  a  b  mm 


(A-l) 


where 


t 

K 


SF 


N 

m 


E 

m 


effective  received  power,  watts 

2 

target  element  radiance,  watts/m  -ster. 

2 

ft(Tr/4)  d  p  (a  "receiver  factor") 

effective  average  background  radiance  contributing  to  scattering, 
watts/m^-ster . 

scattered  fraction;  the  fraction  of  "average"  background 

vR 

radiance  scattered  into  the  radiometer  FOV.  SF  =  (o/y(l  -  e  ) 

blackbody  radiance  at  rhe  air  temperature  of  the  transmission  path 

the  effective  emissivity  of  the  air  mass  within  the  radiometer 
FOV.  The  value  of  e  is 


=  detector  element  FOV,  steradians 
d  =**i  receiver  effective  clear  aperture  diameter,  meters 
n  =  detector  quantum  efficiency 

a  =  atmospheric  scattering  coefficient,  per  kilometer 

k  =  atmospheric  absorption  coefficient,  per  kilometer 

Y  =  atmospheric  attenuation  coefficient,  per  kilometer 

T  =  path  transmittance,  =  e 

4 
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PATH  EQUATION  NORMALIZED  TO  RADIANCES 


Equation  (A-l)  shows  the  received  power  to  be  composed  of  a  target 
component,  an  atmospheric  scattering  component,  and  an  air  mass  radiance 
component,  respectively.  If  the  receiver  factor,  R,  is  factored  out, 
Equation  (A-l)  may  be  written  as: 


N 

a 


NT  +  N,  (SF)  +  N  e  * 
tab  mm 


(A- 2) 


where  N^  is  an  "apparent"  or  measured  radiance.  Hence  the  measured  radiance 
is  composed  of  target,  atmospheric  scattering,  and  atmospheric  emission 
radiance  components. 


SCOPE  OF  APPLICATION 

Equations  (A-l)  and  (A-2)  apply  strictly  only  to  a  uniform  horizontal 
path  and  to  narrow  spectral  regions  over  which  the  absorption  and  scattering 
coefficients  are  constant,  so  that  exponential  transmittance  applies.  Under 
other  circumstances,  more  accurate  estimates  may  be  obtained  through  use  of 
one  of  the  available  band  transmittance  models  (e.g.,  LOWTRAN) .  Such  models 
use  experimentally-derived  empirical  relationships  for  calculating  the 
average  optical  absorption  of  atmospheric  constituents  which  have  highly 
structured  spectral  characteristics. 


With  appropriate  caution,  and  with  due  regard  for  the  approximate 
nature  of  the  results,  Equations  (A-l)  and  (A-2)  have  been  found  useful  in 
estimating  horizontal  path  transmittances  and  radiances  and  in  pointing  the 
way  to  calculation  of  slant  path  transmittances  and  radiances  for  wavelength 
bands. 

VERTICAL  AND  SLANT  PATH  CALCULATIONS 

For  vertical  and  slant  path  calculations,  the  first  term  of  Equation 
(A-2),  NT  ,  is  relatively  easy  to  calculate.  The  scattering  and  thermal 

t  cl 

*  This  expression  is  a  much-simplified  version  of  La  Rocca's  general  equation 
for  "the  atmospheric  radiative  transfer  equation,"  Equation  (1)  of  Ref.  17. 
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radiance  terms,  however,  must  ideally  be  calculated  on  a  continuously  variable 
or  layer-by-layer  basis,  with  the  radiance  of  each  differential  thickness  or 
layer  attenuated  to  the  observer  location  by  the  transmittance  of  the  inter¬ 
vening  path.  Such  computations  can  be  tedious  and  time-consuming,  even  when 
layers  of  one-kilometer  thickness  are  assumed,  as  in  using  the  McClatchey 
coefficient  tables  from  References  1,  4,  and  5. 

It  was  determined,  however,  that  the  integrated  scattering  fraction, 

SF,  and  the  integrated  air  mass  emissivity,  e  ,  over  a  vertical  or  slant  path 

m 

are  independent  of  the  variation  of  a  and  k  with  altitude.  Hence,  if 

and  were  reasonably  constant  with  altitude,  the  layer-by-layer  computa¬ 
tion  would  be  unnecessary.  This  simplification  is  reasonable  for  N,  ,  but  not 

b 

for  N^,  hence  its  applicability  is  essentially  limited  to  the  visible  band. 

DERIVATION  OF  SIMPLE  PATH  RADIANCE  FORMULA  FOR  VISIBLE  WAVELENGTHS 

Taking  the  layer-by-layer  approach,  the  path  radiance,  looking 
vertically  downward  from  a  point  above  the  earth,  is  given  by: 


£ 

all  Ar 


T  ff/Ar  \  -Y(Ar  )\ 

N  i)  1  -  e 

b(Ar  )  Y(Ari) 

i  L  /  , 


-ET(r.) 


(Ar^ 

"<4ri) 


-y(Ari)\  /  -EyC^) 


(A-3) 


where  the  first  term  is  the  scattering  component  and  the  second  is  the  thermal 

til 

radiation  component.  The  (Ar^)  subscript  signifies  the  i  layer  of  thick¬ 
ness  Ar,  and  Ey.^^  signifies  the  total  optical  depth  from  the  observer  to 
the  i^  layer.  For  a  slant  path,  each  a,  k,  and  y  is  multipled  by  sec  0, 
where  0  =  path  zenith  angle. 


In  the  visible  wavelength  region,  absorption  is  negligible  in  comparison 
to  scattering,  while  in  the  8-13  micrometer  IR  region,  scattering  is  small 


compared  to  absorption  and  resultant  thermal  emission.  Hence  Equation  (A-3) 
can  be  considerably  simplified  for  computation  in  either  of  these  regions. 

In  the  visible  region,  for  example,  it  can  be  reduced  to 


N  Z  [  N 

all  Ar.  (Ar  ) 

1  L  i 


Then,  if  N,  is  considered  constant  with  altitude,  Equation  (A-4)  for  the 
b 

visible  region  reduces  to 
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all  Ar , 


■ 

,  -a (Ar . ) 

-Za(r  ) 

1  -  e  l 

e  i 

• 

(A-5) 


f 


i 


Dividing  through  by  and  writing  the  summation  term  by  term,  Equation 

(A-5)  becomes 


Scattered  fraction 
contributions  by 
layers,  assuming 
each  layer  is 
1  km  thick  and 
a  values  are  in 
km~l.  Also,  that 
observer  is  12  km 
above  surface. 


Performing  the  indicated  multiplications: 


-o  -(a  + 

SF  =  1  -  e  A  e 


•  -o12) 


■  -  I 


-  (<?3  + 


'a12) 


-°3  +  ’  •  -°12> 
+  1  -  e  e 
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- (°2  +  ■  ■  •  «12)  -(»!  +  »2  +  •  '  '  a12) 

SF  =  e  -  e 


+e 


(<,3  +  '  '  '  °12>  -(o2  +  «3  +  •  •  '  "12> 


-  e 


(aA  +  *  •  *  a19)  -(oQ  +  aA  +  •  ■  •  a12) 


+e 


12  3  4 

-  e 


+e 


12 


-  e 


-<0U  +  "121 


+  1 


-  e 


12 


1  -  e 


-(ax  +  a2  + 


a12> 


=  1  -  T 


This  Is  the  simple  formula  for  path  radiance  (scattered  fraction)  that  is 
commonly  used  for  horizontal  path  calculations  at  visible  wavelengths.  Switch¬ 
ing  to  photometric  terminology,  the  visible  path  brightness,  B  ,  is  given  by: 

s 

B  =  B  (1  -  T  )  for  any  path  over  which  the  average  background 

S  D  cl 

brightness,  B^,  can  be  assumed  reasonably  constant. 

A  check  on  this  result  was  performed  using  the  layer-by-layer 
approach.  The  results  are  tabulated  in  Table  A-l  (for  the  SIS  66-km 
slant  path,  and  a  mid-latitude/summer  "clear"  atmosphere).  The  reader  is 
also  referred  to  Bolleau  (Ref.  18)  for  tables  of  slant  path  radiances 
(luminances),  measured  from  an  aircraft  in  flight. 


IR  WAVELENGTH  REGION  OF  8-13  MICROMETERS 

In  this  wavelength  region,  scattering  is  small  compared  to  absorption 
and  resultant  thermal  emission,  so  that  Equation  (A-3)  can  be  reduced  to: 


z 

all  Ar 

m 

Nm(Ari) 

(l  -  e'k(^)l 

e'Ek<ri>  |  ' 

i  /  - 

1 

- 

1  / 

(A-6) 
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Table  A-l.  Example  Layer  Summation  of  Path  Radiances  at  Visible 
Wavelengths;  Clear  Atmosphere,  66-km  Slant  Path 


1 

2 

3 

4 

5 

Layer 

Ave  o  * 

SF  =  ** 

£  a 

T  to  Layer 

SF 

Layer 

(km) 

(km  1) 

~aat 

1  -  e 

t 

to  Layer 

-Ea  a 

=e 

Contrib . 

=  3x5 

0-1 

0.1263 

0.495 

0.19893 

0.341 

0.169 

1-2 

0.0615 

0.283 

0.13743 

0.475 

0.134 

2-3 

0.0323 

0.160 

0.10513 

0.566 

0.091 

3-4 

0.02026 

0.104 

0.08487 

0.632 

0.066 

4-5 

0.01567 

0.0813 

0.06920 

0.688 

0.056 

5-6 

0.01305 

0.0682 

0.05615 

0.738 

0.050 

CO 

CM 

uO 

6-7 

0.01132 

0.0594 

0.04483 

0.785 

0.047  ' 

CM 

CO 

7-8 

0.01045 

0.0550 

0.03438 

0.830 

0.046 

o  * 
* 

8-9 

0.00970 

0.0511 

0.02468 

0.875 

0.045 

CM 

• 

9-10 

0.00892 

0.0471 

0.01576 

0.918 

0.043 

lO  rH 

1  a)  a 

10-11 

0.00817 

0.0432 

0.00759 

0.960 

0.041 

ii  H 

11-12 

0.00759 

0.0402 

0 

1 

0.040 

(T3 

H 

Totals 

0.32523 

0.827*** 

* :  a  = 

t 

ct  +  c  at 
m  a 

0.5145  micrometer,  from 

AFCRL-72-0497 . 

**:  a  =  sec  6  =  5.41 


***:  The  resulting  total  SF  =  0.827,  which  agrees  with  1  -  T  ,  =  1  -  0.172 
=  0.828.  3 


0.172*** 


^m(,'r  )  t^*e  e1u^va^ent:  blackbody  radiance  of  the  i  Ar  layer,  and  its 

value  is  a  moderately  strong  function  of  altitude.  Hence  the  simplification 

made  at  this  point  for  visible  wavelengths  (N,  -  constant  with  altitude)  has 

D 

no  counterpart  in  this  wavelength  region,  and  the  path  radiance  is  not 
approximately  given  by  a  constant  multiplied  by  (1  -  T  ) . 

Si 

Layer-by-layer  summation  of  path  thermal  radiance  (for  the  mid¬ 
latitude  summer  model)  was  calculated  at  10.591  micrometers  for  the  SIS 
66-kilometer  slant  path,  as  shown  by  Table  A-2.  The  results  were 

2 

Nt  =  1 5 . 9/ 7t  =  5.06  watts/m  -ster-pm,  looking  down. 


and 


2 

=  27.0/ir  =  8.59  watts/m  -ster-pm,  looking  up. 

It  is  clear  from  this  example  that  the  path  radiance  (thermal)  can 

not  be  given  by  N  (1  -  T  )  since  the  path  radiance  is  different  "looking  up" 
m  3 

and  "looking  down,"  even  though  T  must  be  the  same  in  both  directions. 

Si 

The  layer-by-layer  summation  process  was  exercised  to  obtain 
approximate  average  values  of  path  radiance  and  path  transmittance,  for  the 
8-12  micrometer  band.  Hite  process  was  repeated,  for  "clear"  and  "hazy" 
atmospheres,  for  two  latitude/season  models,  and  for  5  horizontal  ranges 
(0,  9,  18,  27,  and  35  nautical  miles),  looking  down  from  40,000  feet  altitude. 
The  starting  points  for  these  calculations  were  estimated  average  values  of 
66-km  slant  path  transmittances  read  from  the  computer-calculated  graphs* 

(see  Section  A-6). 

*Note:  L0WTRAN  3  runs  of  August-September  1975  were  used. 

Results  as  shown  in  Table  A-3  and  Figure  A-9  would  be 
modified  somewhat  if  the  computations  were  repeated  using 
"Modified  L0WTRAN  3"  results.  With  respect  to  Figure  A-9, 
a  brief  "eyeball"  comparison  of  L0WTRAN  3  and  Modified 
L0WTRAN  3  indicates  approximate  improvements  a9  follows, 
for  the  35  nm  and  0  nm  (vertical)  transmittances  of 
Figure  A-9  using  modified  LOWTRAN  3;  (note  continued  on 
page  22) 
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Table  A-2.  66-km  Slant  Path  Air  Mass  Thermal  Radiance  at  Receiver 

at  10.591  ym,  Assuming  Molecular  Absorption  Only 


Total  vertical  0D  «  0.94359;  T  =  e  5,41 94359)  =  0>0061 

ci 


Layer 

(km) 

°K 
Ave . 
Temp. 

@  10.6 

N  * 
bb 

W/m^-ym 

Ave 

k 

m 

(knT-*-) 

Air  Mass 
Emiss. 
Layer  e 
-5.41k 

m 

=l-e 

I0D 

to  Layer 

T 

a 

to  Layer 

-5.4170D 

=e 

W/m^-ym 

N  * 

s 

Contrib . 

Looking  down:  mic 

-latitude 

summer 

(very  clear) 

0-1 

292 

27.8 

0. 3256 

0.828 

0.61799 

0.035 

0.806 

1-2 

287.5 

25.7 

0.1877 

0.638 

0.43029 

0.098 

1.60 

2-3 

282 

23.7 

0.1152 

0.464 

0.31509 

0.182 

2.00 

3-4 

276 

21. 3 

0.07582 

0.336 

0.23927 

0.274 

1.96 

4-5 

270 

19.1 

0.05544 

0.259 

0.18383 

0.370 

1.83 

5-6 

264 

16.9 

0.04468 

0.215 

0.13915 

0.471 

1.71 

6-7 

258 

15.0 

0.03752 

0.184 

0.10163 

0.577 

1.59 

7-8 

251.5 

13.3 

0.03018 

0.151 

0.07145 

0.679 

1.36 

8-9 

245 

11.7 

0.02378 

0.121 

0.04767 

0.773 

1.09 

9-10 

238.5 

9.72 

0.01952 

0.100 

0.02815 

0.859 

0.83 

10-11 

232 

8.40 

0.01574 

0.082 

0.01241 

0.935 

0.64 

11-12 

225.5 

7.18 

0.01241 

0.065 

0 

1 

0.47 

Total  = 

0.94359 

Total  = 

15.9* 

Looking  up:  mid-latitude  summer  (very  clear) 


0-1 

27.8 

0.828 

0 

1 

23.02 

1-2 

25.7 

0.638 

0.3256 

0.172 

2.82 

2-3 

23.7 

0.464 

0.5133 

0.0622 

0.68 

3-4 

21.3 

0.336 

0.6285 

0.0334 

0.24 

4-5 

19.1 

0.259 

0.70432 

0.0221 

0.11 

5-6 

16.9 

0.215 

0.75976 

0.0164 

0.06 

6-7 

15.0 

0.184 

0.80444 

0.0129 

0.04 

7-8 

13.3 

0.151 

0.84196 

0.0105 

0.02 

8-9 

11.7 

0.121 

0.87214 

0.00893 

0.01 

9-10 

9.72 

0.100 

0.89592 

0.00785 

0.008 

10-11 

8.4 

0.082 

0.91544 

0.00707 

0.005 

11-12 

7.18 

0.065 

0.93118 

0.00649 

0.003 

Total  =  27.0* 


*  Radiant  emittances ; divide  Ng  by  it  to  get  radiances. 
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Mid-latitude  summer  Mid-latitude  winter 

35  run  vertical  35  nm  vertical 

70-80%  10-12%  20-25%  3-5% 

Correction  of  the  path  radiance  curves  has  not  been 
attempted;  however,  some  decrease  in  values  would 
result. 

Estimated  average  values  of  total  path  and  kilometer  layer  values  of 
y,  o,  and  k  were  derived  from  these  computer-calculated  transmittance  values 
and  previously  calculated  values  of  aerosol  attenuation  for  this  path.  It  is 
cautioned  that  this  approach  is  an  approximation,  for  the  reason  previously 
stated.  However,  experience  has  shown  that  the  approximation  is  a  reasonable 
one  for  the  circumstances  of  this  particular  problem.  (The  alternate  approach 
is  a  host  of  LOWTRAN  runs  to  produce  the  layer-by-layer  radiance  and  trans¬ 
mittance  data  required,  or  a  new  path  radiance  computer  program  which  would 
compute  path  radiance  more  directly.) 

Table  A-3  is  a  sample  tabulation  (1  of  16)  involved  in  the  manual 
calculations.  The  model  atmosphere  air  temperatures  are  from  McClatchey 
(Ref.  1),  and  the  associated  blackbody  radiance  values  were  derived  using  the 
G.E.  Radiation  Calculator.  Layer-by-layer  average  molecular  absorption 
coefficient  values,  k^,  were  taken  as  0.303  of  McClatchey's  values  at 
10.591  micrometers,  for  mid-latitude  summer,  and  0.425  for  mid-latitude 
winter,  based  upon  the  cited  derivations,  and  altitude  variations  in  the 
aerosol  absorption  coefficient  were  taken  at  the  same  ratio  as  in  McClatchey. 

IR  WAVELENGTH  REGION  OF  3-5  MICROMETERS 

In  this  wavelength  region  the  atmosphere  exhibits  optical  properties 
which  are  a  composite  of  those  for  the  shorter  and  longer  wavelength  bands. 
Molecular  and  aerosol  absorption  and  aerosol  scattering  are  all  significant 
in  the  attenuation  equation,  and  both  scattering  and  thermal  radiation 
should  be  considered  in  estimating  path  radiance.  Portions  of  this  band 
(in  the  vicinity  of  2.2  and  3.8  micrometers)  exhibit  optimum  optical  trans¬ 
mission  over  a  wide  range  of  atmospheric  variables.  Optical  energy  in  these 
wavelength  regions  penetrates  atmospheric  haze  much  better  than  shorter 
(e.g.,  visible)  wavelengths,  and  is  much  less  affected  by  atmospheric  water 
vapor  than  longer  wavelengths  (e.g.,  the  8-13  micrometer  IR  region). 
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Table  A-3.  Estimated  66-km  Slant  Path  Air  Mass  Thermal  Radiance 
over  8-12  pm  Wavelength  Region  Looking  Down 
from  40  K'  (12.192  km) 

Horizontal  range  35  nm,  64.82  km;  9  =  79.35°;  sec  0  =  5.41  (=a) 


Layer 

T 

oK 

N  * 
bb 

Emissivity 

to  Layer 

W/m^-pm 

N  * 

Layer 

Ave . 

Watts/ 

k 

f- 

“aYt 

0D 

-aEOD 

s 

(km) 

Temp. 

m^-ster-pm 

1  ■ 

Yt 

i1_G 

to  Layer 

e 

Contrib 

_ 

Mid-latitude  summer  (hazy) 


0-1 

292 

8.83 

0.18454 

0.544 

0.29119 

0.207 

0.994 

1-2 

287.5 

8.03 

0.0Q668 

0.  368 

0.19451 

0.349 

1.03 

2-3 

282 

7.42 

0.05515 

0.242 

0.13936 

0.471 

0.846 

3-4 

276 

6.67 

0.03449 

0.164 

0.10487 

0.567 

0.620 

4-5 

270 

5.93 

0.02439 

0.121 

0.08048 

0.647 

0.464 

5-6 

264 

5.10 

0.01941 

0.098 

0.06107 

0.719 

0.359 

6-7 

258 

4.57 

0.01629 

0.083 

0.04478 

0.785 

0.298 

7-8 

251.5 

4.00 

0.01316 

0.068 

0.03162 

0.843 

0.229 

8-9 

245 

3.42 

0.01044 

0.054 

0.02118 

0.892 

0.165 

9-10 

238.5 

2.86 

0.00862 

0.045 

0.01256 

0.934 

0.120 

10-11 

232 

2.49 

0.00699 

0.036 

0.00557 

0.970 

0.087 

11-12 

225.5 

2.09 

0.00557 

0.029 

0 

1 

0.061 

Totals  = 

0.47573 

(T  =  0.076) 
a 

5.27 

Mid-latitude  winter  (hazy) 

0-1 

8.83 

0.06832 

0.193 

0.116415 

0.533 

0.908 

1-2 

8.03 

0.03477 

0.125 

0.081645 

0.643 

0.645 

2-3 

7.42 

0.02107 

0.090 

0.060575 

0.721 

0.481 

3-4 

6.67 

0.01452 

0.069 

0.046055 

0.779 

0.359 

4-5 

5.93 

0.01064 

0.054 

0.035415 

0.826 

0.265 

5-6 

5.10  . 

0.008363 

0.043 

0.027052 

0.864 

0.189 

6-7 

4.57 

0.006842 

0.035 

0.02021 

0.896 

0.143 

7-8 

4.00 

0.005569 

0.029 

0.014641 

0.924 

0.107 

8-9 

3.42 

0.004587 

0.024 

0.010054 

0.947 

0.078 

9-10 

2.86 

0.003600 

0.018 

0.006454 

0.966 

0.050 

10-11 

2.49 

0.003277 

0.017 

0.003177 

0.983 

0.042 

11-12 

2.09 

0.003177 

0.016 

0 

1 

0.033 

Totals  = 

0.184735 

(T  =  0.368) 
a 

3.30 

Columns  4  and  5  data  are  from  a  separate  tabulation. 

*  Radiant  emittances;  divide  N  by  it  to  get  radiances;  N,  *  (Layer  e)  (T  to 

i  \  S  D  D  3 

layer) . 

Note:  The  note  on  page  20  regarding  LOWTRAN  3  applies  to  this  page  also. 
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The  path  radiance  in  this  region  will  be  relatively  small  for  two 
reasons:  (1)  the  natural  daylight  illumination  is  much  reduced  compared  to 

the  visible,  and  (2)  thermal  radiance  is  much  reduced  compared  to  the  8-13 
micrometer  IR  band.  However,  accurate  calculations  of  path  radiance  are  much 
more  laborious  because  the  simplifying  approximations  of  those  bands  do  not 
apply  in  the  3-5  micrometer  region. 

We  have  not  attempted  to  make  estimates  of  path  radiance  for  the 
66-kilometer  SIS  slant  path  in  this  wavelength  region.  However,  computer 
calculated  band  transmittances  and  manually  calculated  transmittance  values 
for  many  DF  laser  lines  in  this  region  have  been  accomplished,  as  reported 
in  Sections  A-6  and  A-10. 

AEROSOL  MODELS 

The  aerosol  model  used  in  the  calculations  reported  herein  is  AFCRL's 
"average  continental  aerosol  model"  (1974)  supplied  to  SRL  by  Reference -19* 
(More  recently,  this  model  appears  in  the  aerosol  spectral  data  table  of  the 
LOWTRAN  3  computer  code,  Reference  8,  p.  86.  The  LOWTRAN  3  report,  dated  7 
May  1975,  was  distributed  in  December  1975.)  The  vertical  distribution  of 
particle  densities  for  this  model  is  given  on  page  9  of  Reference  1,  for 
normalized  "clear"  (S.L.  visibility  =  23  km)  and  normalized  "hazy"  (S.L. 
visibility  =  5  km)  atmospheres.  In  SRL’s  normalized  "light  haze"  atmosphere 
(S.L.  visibility  =  10.8  km)  the  particle  density  at  each  altitude  is  the 
geometric  mean  of  the  "clear"  and  "hazy”  values. 

AFCRL  also  has  an  "estimated  marine  aerosol  model,"  an  advance  copy 
of  which  was  supplied  to  SRL  by  Reference  20,  and  several  other  models 
(i.e.,  urban  and  rural)  intended  for  near-future  publication  in  a  supplement 
to  LOWTRAN  3.  We  have  exercised  the  marine  model  to  some  extent  on  the 
SIS  study,  as  reported  in  other  sections  of  this  appendix. 

RESULTS 

Some  typical  calculated  results  are  illustrated  in  Figures  A-5 
through  A-9,  and  others  may  be  found  elsewhere  in  this  appendix. 
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Figures  A— 5  through  A- 7  show  calculated  values  of  slant  path  trans¬ 
mittance  at  wavelengths  of  0.6525  micrometer,  0.55  micrometer,  3.73  micro¬ 
meter,  and  10.591  micrometer.  The  visible  wavelength  curves  are  nearly 
independent  of  latitude-season  model.  For  example,  the  curves  for  0.55 
micrometer  and  0.6525  micrometer  apply  within  0.002  for  any  of  the  five 
McClatchey  latitude-season  models.  On  the  other  hand,  the  large  dependence 
upon  haze  level  is  clearly  shown,  for  these  wavelengths. 


At  the  other  extreme,  the  small  dependence  of  10.591  micrometer 
transmittance  upon  haze  level,  and  the  large  dependence  upon  latitude-season 
model,  are  clearly  shown.  The  3.73  micrometer  curves  show  an  intermediate 
level  of  dependence  upon  both  factors. 

Figure  A-8  illustrates  slant  path  luminance  at  0.55  micrometer,  the 
mid-point  of  the  visible  band,  calculated  from  the  simple  formula  previously 
discussed. 

Figure  A-9  illustrates  calculated  band  transmittance  averages  and 
path  radiances  for  the  8-12  micrometer  band,  computed  according  to  the 
layer-by-layer  method  previously  described.  Both  sets  of  curves  (trans¬ 
mittance  and  path  radiance)  show  the  cited  strong  dependence  upon  latitude- 
season  model,  with  lesser  dependence  upon  the  haze  level.  Comparison  of 
these  curves  with  those  of  Figure  A-7  shows  the  relatively  better  trans¬ 
mittance  of  the  8-12  micrometer  band  as  compared  to  monochromatic  trans¬ 
mittance  of  a  CC^  laser  wavelength.  (Also,  see  Note,  page  20.) 

While  comparable  curves  comparing  3-5  micrometer  band  transmittance 
with  DF  laser  monochromatic  transmittances  have  not  been  developed,  review 
of  available  data  shows  that  most  DF  laser  lines  exhibit  better  atmospheric 
transmittance  than  the  average  band  transmittance.  This  is  opposite  to  what 
occurs  in  the  8-12  micrometer  band  with  respect  to  band  and  CO2  laser  line 
transmittances . 
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Hazy 


-6.  Typical  Calculated  Slant  Path  Transmit tances 
iQO  Feet  to  Sea  Level,  3.73  and  0.55  Micrometer 


Figure  A-7.  Typical  Calculated  Slant  Path  Transmit tances 
40,000  Feet  to  Sea  Level,  10.591  pm 


Typical  Calculated  Path  Luminance  for  40,000  Feet  to  Sea  Level  Path 
Sky  Brightness  =  3400  Lumens/m  -ster ,  Visible  Wavelength 
Region  (0.55  pm),  Mid-Latitude  Summer 
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Section  A-4 

ATMOSPHERIC  OPTICAL  TURBULENCE 


INTRODUCTION 

Identification  schemes  which  utilize  optical  signals  will  be  adversely 
affected  by  the  propagation  of  these  signals  through  a  turbulent  atmosphere. 

It  is  convenient  to  divide  the  propagation  path  into  two  parts,  the  turbu¬ 
lent  aircraft  boundary  layer  and  the  free  atmosphere.  The  reason  for  this 
division  is  that  the  free  atmosphere  optical  turbulence  is  caused  primarily 
by  thermally  induced  air  density  fluctuations,  and  is  characterized  by  a 
Kolmogorov- type  refractive  index  fluctuation  spectrum,  while  the  boundary 
layer  optical  turbulence  is  caused  primarily  by  mechanically  induced  density 
fluctuations  whose  structure  is  quite  different. 

The  theory  of  free  atmosphere  optical  turbulence  has  been  developed 
and  experimentally  verified  sufficiently  to  permit  the  direct  calculation 
of  optical  signal  degradation.  However,  the  theory  of  aircraft  boundary 
layer  optical  turbulence  is  not  so  well  developed,  and  experimental  results 
are  only  presently  being  obtained.  For  the  present  purpose,  we  must  rely 
on  preliminary  data  for  estimates  of  the  aircraft  boundary  layer  effect. 

Since  the  types  of  degradations  produced  by  the  boundary  layer  are  the 
same  as  those  produced  by  the  free  atmosphere,  we  proceed  to  a  brief  des¬ 
cription  of  free  atmosphere  turbulence. 

FREE  ATMOSPHERE  TURBULENCE 
Background 

The  origin  of  the  free  atmosphere  turbulence  effects  lies  in  the 
fact  that  the  refractive  index  at  a  point  in  the  atmosphere  exhibits  small 
fluctuations  about  its  mean  value.  The  fluctuations  are  correlated  only 
over  short  distances,  and  hence  produce  a  general  degradation  in  optical 
signals  propagating  through  the  atmosphere.  Clearly,  a  detailed  descrip¬ 
tion  of  the  degradation  is  not  possible  without  a  detailed  knowledge  of  the 
refractive  index  space-time  dependence.  However,  a  statistical  description 
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of  the  refractive  index  fluctuations  does  permit  a  statistical  description 
of  the  resultant  optical  degradation.  It  is  this  approach  which  has  been 
used  by  many  authors  to  obtain  useful  formulae  for  calculating  the  various 
degradations  produced  by  atmospheric  turbulence.  A  major  portion  of  the 
theory  is  contained  in  the  volume  by  Tatarski  (Ref. 14),  and  the  review 
articles  by  Lutomirski  et  al(Ref,13)  and  by  Lawrence  and  Strohbehn  (Ref .21) 
provide  useful  practical  discussions  and  summaries  of  more  recent  important 
theoretical  and  experimental  results. 

Refractive  index  fluctuations  arise  primarily  from  temperature 

fluctuations  in  the  atmosphere,  which  in  turn  are  produced  by  the  turbulent 

mixing  of  hot  air  (usually  heated  by  the  earth's  surface)  with  cool  air. 

The  mixing  occurs  through  a  series  of  eddies  ranging  in  size  from  Lq 

(-  1  meter  near  the  surface)  to  Z  (-1  millimeter).  L  and  S,  are  called 

o  o  o 

the  outer  and  inner  scales  of  turbulence.  Pressure  and  humidity  fluctuations 
also  produce  index  fluctuations,  but  their  effect  is  usually  negligible 
compared  with  temperature  fluctuations.  An  exception  to  this  rule  occurs 
when  dry  air  mixes  with  moist  air,  as  might  occur  at  a  land-sea 
interface  (Ref. 22). 

A  useful  descriptor  of  the  index  fluctuations  at  a  point  z  along  a 
propagation  path  is  the  spatial  index  power  spectrum  ♦  ^(*£,z),  which  is  defined 
as  the  three-dimensional  Fourier  transform  of  the  index  spatial  correlation 
function.  For  locally  isotropic  turbulence,  $n(£,z)  is  often  written 
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(A- 7) 


I 


where  K  is  the  spatial  frequency  and  (z)  is  a  parameter  which  is  a 

measure  of  the  strength  of  the  turbulence  (the  rms  index  fluctuation  is 
1/3 

approximately  )  •  The  gross  spatial  dependence  of  is  contained 

in  C^Cc).  For  many  applications,  the  dependence  on  SL^  and  Lq 

can  bo  supressed,  yielding  the  Kolmogorov  spectrum 

?  -11/3 

<5>  (K, z)  =  0.033  C  (z)  K  (A-8) 

n  n 

Since  no  good  models  exist  for  the  outer  scale  at  high  altitudes  and  since 

the  inner  scale  can  usually  be  neglected,  we  utilize  the  Kolmogorov  spectrum 

2  2/3 

for  this  application.  We  note  that  C  has  the  dimension  of  (length) 

n 

and  that  the  constant  0.033  is  appropriate  when  MKS  units  are  used. 

A  useful  conceptualization  of  the  refractive  index  field  is  to 
imagine  a  dynamic  collection  of  arbitrarily  shaped  lenses  having  different 
refractive  indices  with  sizes  ranging  from  to  L^.  An  optical  wave 
passing  through  the  atmosphere  will  then  undergo  a  series  of  focusings 
and  steerings,  and  will  impinge  on  a  target  or  receiver  with  considerable 
degradation.  The  degradation  will  be  manifest  in  a  variety  of  ways,  de¬ 
pending  on  the  type  of  optical  wave  involved  (i.e.,  spherical  wave,  focused 
beam,  or  reflected  ambient  light),  the  viewing  techniques  (i.e.,  point 
detection  or  imaging  system) ,  and  the  coherence  of  the  initial  optical 
signal  as  well  as  other  factors. 

Most  treatments  divide  the  total  wave  degradation  for  infinite 
extent  wavefronts  into  two  parts:  amplitude  fluctuations  and  phase  fluc¬ 
tuations.  The  amplitude  fluctuations,  termed  scintillation,  are  produced 
primarily  by  the  focusing  action  of  the  index  field  and  are  observed  as  a 
breakup  of  the  wavefront  into  fluctuation  patches  of  high  and  low  intensity. 
The  phase  fluctuations  are  produced  primarily  by  the  steering  action  of 
the  index  field  and  are  observed  to  first  order  as  fluctuations  in  the 
angle  of  arrival  (wavefront  propagation  direction)  of  the  radiation. 

e 

I 


The  same  division  (amplitude  and  phase)  is  generally  carried  over 
in  treatments  of  finite  beam  degradation.  However  the  turbulence  will  also 
produce  two  additional  effects,  motion  of  the  beam  centroid  (beam  wander) 
and  fluctuations  and  enlargement  of  the  beam  waist  (beam  spread).  These 
effects  are  the  result  of  first  order  and  higher  phase  fluctuations 
respectively.  In  passing,  we  note  that  intensity  fluctuations  actually 
observed  at  a  point  in  the  beam  may  be  due  both  to  scintillation  and  to 
beam  wander.  The  scintillation  which  we  estimate  below  should  be  taken  as 
the  intensity  fluctuation  at  the  beam  center. 

In  imaging  systems,  the  wavefront  incident  on  the  receiver  will  be 
degraded  by  scintillation  and  by  phase  distortion.  However,  it  is  not 
these  degradations  directly  but  rather  those  observed  in  the  image  plane 
which  are  of  interest.  In  general,  a  point  in  the  image  will  show  power 
fluctuations  due  to  scintillation,  motion  (image  dancing)  due  to  first 
order  phase  fluctuations,  and  spreading  (increase  in  the  point  spread 
function  (PSF)  width)  due  to  higher  order  phase  degradations.  The  power 
fluctuations  are  often  negligible  due  to  the  averaging  effect  of  the  re¬ 
ceiver  and  source. 

Before  discussing  the  specific  turbulence-produced  effects  studied 
here,  we  point  out  that  the  time  frame  in  which  observations  are  made  can 
be  important.  For  example,  in  a  short  time  (approximately  10  ms)  imaging 
system,  the  resolution  will  be  determined  by  the  instantaneous  point  spread 
function,  and  will  not  be  affected  by  the  image  dancing.  However,  for  long¬ 
time  averaging  (as  in  a  10  second  photographic  exposure),  the  image  dancing 
will  be  convoluted  with  the  PSF,  yielding  a  further  decrease  in  resolution. 
When  each  image  point  is  a  short-time  average,  but  the  image  frame  is 
collected  over  a  long  time,  the  overall  effect  will  be  a  high-resolution 
image  (as  determined  by  the  instantaneous  PSF) ,  with  considerable  distortion 
as  determined  by  the  image  dancing. 
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Analysis 

A  general  optical  ship  identification  system  will  potentially 
utilize  an  illuminator  and  an  imaging  receiver.  The  turbulence  degradations 
which  need  be  considered  for  the  illuminator  are  scintillation,  beam  wander, 
instantaneous  beam  spread,  and  long-term  beam  spread.  For  an  imaging 
receiver,  we  need  to  consider  scintillation,  image  dancing,  instantaneous 
point  spread  function,  and  long-term  point  spread  function.  We  now  present 
the  formulae  used  to  estimate  these  effects.  Since  formulae  for  finite 
beam  propagation  are  complex  and  unwieldy,  and  since  plane  wave  or  point 
source  (spherical  wave)  formulae  often  yield  adequate  approximations,  we 
use  those  formulae  where  applicable. 

2 

The  descriptor  for  scintillation  is  a  ,  the  variance  of  the  dis- 

X* 

tribution  for  the  natural  logarithm  of  the  normalized  intensity  as  observed 

by  a  point  receiver.  For  spherical  wavefronts  of  wave  number  k  =  2tt/X, 

2 

propagating  from  position  z'  =  0  to  position  z'=z,o  is  given  by 
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For  the  Kolmogorov  spectrum,  this  equation  reduces  to 
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The  corresponding  equation  for  plane  waves  is 


a.  =  2.24  k7/6  /  (z-z')5/6  C  2(z’)dz'  (A-ll) 

^’P  o  n 

We  note  that,  for  spherical  waves,  turbulence  near  the  path  midpoint  is 
weighted  strongest,  but  for  plane  waves  the  strongest  weighting  is  for 
turbulence  near  the  source. 

2 

These  equations  yield  good  agreement  with  measurements  for  a <  1 ; 

however,  for  strong  turbulence  and/or  long  propagation  paths  such  that 
2 

>  1,  the  formulae  yield  values  which  are  too  large.  Experimentally,  a 

saturation  of  scintillation  occurs.  An  empirical  correction  factor  has  been 

obtained  (Ref. 23)  which  yields  a  corrected  a  given  by 

x, ,  c 


a 


i,  c 


1  +  Ao* 


(A-12) 


where  A  and  B  are  constants  with  a  slight  wavelength  dependence. 


In  addition  to  the  saturation  correction,  the  averaging  effect  of 
the  receiver  (or  resolution  element  for  the  illuminator)  must  be  considered. 
A  useful  engineering  formula  for  the  aperture  average  correction  factor 
is  (Ref. 13) 


(A-13) 


Then  a. 

Z,av 


where  D  is  the  receiver  diameter  and  is  the  transverse  amplitude 

correlation  length  and  is  given  approximately  by 
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£ 


(A-14) 


This  factor  "a"  is  somewhat  optimistic  for  >1,  but  is  accurate 

enough  for  the  work  at  hand  (Refs  24,  25). 


For  imaging  systems,  one  further  correction  is  needed  to  account  for 
the  averaging  over  the  resolution  element  at  the  target.  Lutomirski  (Ref. 13) 
gives  this  factor  as 


S  = 


—  I  ;  then 

(  2  _  2 
a.  =  S  o„ 

r  / 

|  £  ,s  £  ,c 

(A-15) 


where  is  as  defined  above  and  r  is  the  linear  dimension  of  the 

resolution  element  at  the  target. 


The  primary  descriptor  for  the  phase  degradation  is  the  phase 

structure  function  D  (p )  defined  by 

<P 


D  (p)  =  <[4>(r)  -  4(T*+  t))2> 


(A-16) 


where  i^(r)  is  the  phase  at  point  r  in  the  target  or  receiver  plane  and 
p  is  a  displacement  transverse  to  the  propagation  path.  The  brackets 
indicate  a  time  average.  For  spherical  wavefronts,  D^(p)  is  given  by 


D.  (p,z)  =  8  r2k2  fZ  dz'  /  [1  -  J,,^2-^-)]  *  (K.z')KdK  (n-17) 

<t>*s  0  o  0  z  n 


For  the  Kolmogorov  spectrum,  this  reduces  to 
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D.  a(p>z) 

<P»  S 


=  2.91  k2  p5/3  /  C  2(z')  (— 1  /  dz' 

a  n  l  7  l 


(A-18) 


The  corresponding  formula  for  plane  waves  is 


D.  >>z) 
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2  5/3 
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(A-19) 


Another  useful  descriptor  for  the  wave  degradation  is  the  mutual 
coherence  function  M(p)  which  can  be  defined  as  the  cross  correlation  of 
the  complex  radiation  field  normalized  to  its  vacuum  value  (Ref.  13). 

M(p)  can  be  written 


M(p)  =  e 


-1/2  D  (p) 
<P 


(A-20) 


where  D,(p)  is  the  phase  structure  function. 
<P 


The  separation  pQ  in  the  receiver  plane  at  which  M(po) 
is  called  the  transverse  coherence  length  and  is  given  by 


-1 


D,(p  ,z)  =  2 
<p  o 


(A-21) 


Formulae  (A-20)  and  (A-21)  are  valid  for  all  three  waveforms  (plane,  spheri¬ 
cal  and  beam),  provided  the  appropriate  phase  structure  function  is  used. 

Using  the  Huygens-Fresnel  approach,  Lutomirski  (Ref.  13)  has  shown 
that  the  total  long-term  beam  spread  for  a  Gaussian  beam  can  be  approxi¬ 
mately  written  as 


9  =  (e02  +  eg2) 


1/2 


(A-22) 


where  0q  is  the  diffraction  limit  spread  and  0g  is  caused  by  turbulence 
and  can  be  written  as 
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( A-  2  3 ) 


Here,  p  is  the  transverse  coherence  length,  and  Equation  (A—  23 )  gives  the 
0  -1 

value  appropriate  for  computing  the  e  diameter.  Since  the  Huygens- 
Fresnel  approach  is  used  in  the  derivation  of  Equation  (A-22) ,  the  spherical 
wave  structure  function  is  used  to  compute  p  . 


The  short-term  beam  wander  $  for  homogeneous  turbulence  (C 


constant)  is  given  by  Chiba  (Ref.  26)  to  be 


.2  01  s  r  2  5/6  ul/6 

4>  =  21 . 5  C  z  k 

n 


(A-24) 


-2 

where  $  is  defined  as  the  e  diameter  of  the  beam  wander  distribution. 

This  formula  should  be  approximately  valid  for  inhomogeneous  turbulence  if 

2 

the  appropriate  average  value  of  C  is  used.  Since  beam  wander  is  a  phase 

n 

effect,  we  have  used  the  phase  structure  function  as  a  basis  for  computing 

2 

the  average  ,  obtaining 


z  ,  5/3 

f  c/(z')  f-  dz 

-2  _  8  0  _ _ 

rs 

n  3  z 


(A- 25) 


The  spherical  waveform  of  D,  was  used  because  it  yields  a  more  pessimistic 

-  2  * 
value  for  C 


Reliable  predictive  formulae  for  the  short-term  beam  spread  y  are 

not  available,  so  we  estimate  this  parameter  by  taking  the  root  square 

difference  between  6  and  $ 

s 


2  2  ' 
(6s  -  O 


(A-26) 


This  estimate  is  likely  to  break  down  for  0g  r  3>  due  to  the  different 
approximations  involved  in  obtaining  $  and  0  . 
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The  short-term  rms  image  dancing  is  equal  to  the  rms  angle  of 
arrival  fluctuation  a,  which  is  given  by  (Ref.  13) 


a 


A>T('d) 

A  . . 

kD 


( A—  27) 


where  D  is  the  receiver  diameter.  The  spherical  wave  structure  function 
is  used  here. 

Combining  Equations  (10a)  and  (10b)  of  Ref.  27,  and  applying  a  point 
source  function,  the  long-term  point  spread  function  PSF  is  given  by: 

PSF(x)  =  C  e_0tZ  /  /  M  (£)  M  (p)  e"1(kPx  cos  0)/fpd<J>dp  (A-28) 

0  0  L  u  3 

where  f  and  D  are  the  lens  _ocal  length  and  diameter,  M  (p)  is  the 

s 

spherical  wave  mutual  coherence  function,  and  M^(p/D)  is  the  lens  transfer 
function 


VP/D)  =  |  [cos'1^)  -  (g)  \Jl  -  (§) 


(A- 29) 


-az 


The  factor  e  represents  transmission  losses. 


The  short-term  PSF  width  is  approximated  as  the  root  square 
difference  between  the  long-term  PSF  width  and  the  image  dancing  parameter 
a. 


In  carrying  out  the  detailed  computations,  it  is  necessary  to  adopt 

2 

a  form  for  the  altitude  variation  of  .  No  adequate  theory  exists  to 
predict  a  form  except  at  low  altitudes  and  under  specific  meteorological 
conditions;  we  therefore  rely  heavily  on  measured  profiles.  Figure  A-10 
shows  results  obtained  from  AFVTL  (Ref.  28  ),  summarizing  airplane  measure¬ 
ments  made  on  several  flights.  Figure  A-ll  shows  a  sample  profile  measured 

2 

from  an  ascending  balloon  (Ref.  29).  profiles  have  also  been  inferred 

from  stellar  scintillation  measurements.  In  arriving  at  a  composite  profile, 
many  authors  (Refs  30,  31,  32)  use  terms  such  as 


40 


Typical  profiles  of  re.fractive-i_ndex  structure  coefficient  Cn(h), 
temperature  T,  and  wind  speed  |v|  vs.  altitude  h  above  mean  sea  level. 


V. 


h  (km) 


2 

Figure  A-ll.  Balloon  Cn  Profile  Measurements 
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The  profile  chosen  for  this  effort  is  of  the  form 


(h) 


„  , -2/3  -h/1000  ,  .  ..-17  -2/3 

Kh  e  4-5x10  m 

o 


(A- 31) 


Here,  Kq  represents  the  ground  Jevel  (h  =  1  m)  turbulence  and  the  -2/3 

exponent  is  the  low  level  (h  <  100  m)  altitude  dependence  observed  under 

neutral  conditions  over  water.  (Under  meteorologically  unstable  conditions, 

a  -4/3  exponent  is  more  appropriate;  we  have  chosen  the  more  conservative 

exponent.)  The  value  of  h^  =  1000  m  describes  the  more  rapid  decrease  of 

the  ground  level  influence  at  higher  altitudes.  The  constant  term  5  x 
-17  -2/3 

10  m  is  a  worst-case  estimate  for  the  highly  variable  upper  tropo- 

2  2 
sphere  level.  Ground  level  measurements  generally  find  varying 

from  ~  10~15  m-2/3  or  less  to  5  x  10~13  m-2/30ver  land,  depending  on  the 

local  meteorological  conditions.  Over  sea,  the  upper  limit  is  an  order  of 

magnitude  lower  for  comparable  solar  flux  and  air  temperature  conditions 

■j  r  i  f  1  Q 

(Ref.  33).  We  have  chosen  Kq  values  of  10  ,10  ,  and  10  for  the 

computations . 

Results 

Programs  for  the  Hewlett-Packard  9820  calculator  were  written  to 
evaluate  the  integrals  involved  in  the  formulae.  Calculations  were  per¬ 
formed  for  three  levels  of  turbulence  (K  =  10  10  ^  ,  10  ^)  and  at 

o 

four  wavelengths  (A  =  0.5145,  1.06,  3.80,  and  10.6  micrometer)  covering  the 
ranges  of  primary  interest.  For  the  calculations,  the  transmitter/receiver 
was  assumed  to  be  at  12.2  km  altitude,  with  the  target  at  a  slant  range  of 
65.9  km. 


For  calculating  scintillation  in  the  transmitted  beam,  we  have  used 
the  spherical  wave  approximation.  This  formula  should  be  quite  good  for 
diverging  beams,  and  the  error  should  be  small  for  collimated  beams  for 
the  air-to-ground  path  since  the  strong  turbulence  is  near  the  ground.  (For 
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collimated  beams  propagating  from  ground  level  to  the  aircraft,  the  plane 
wave  approximation  would  be  more  accurate.)  Since  we  are  interested  in 
the  ultirn? te  effect  of  the  scintillation  on  the  target  image,  we  have 
averaged  the  scintillation  over  a  one-meter  spot  corresponding  to  a  15- 

microradian  resolution  element.  Table  A-4  shows  the  results,  including  the 

2  2 
theoretical  a  ,  the  value  corrected  for  saturation  o  (m) ,  and  the  aper- 

2 

ture-averaged  value  a  (A)  ,  along  with  the  amplitude  correlation  length 
and  the  resultant  rms  intensity  fluctuation  level  both  at  a  point  (dB(pt)) 
and  averaged  over  the  one-meter  spot  (dB(avg)). 

Table  A-5  shows  the  short-term  beam  wander  $  for  the  air-to-ground 

2 

path.  Also  shown  is  the  average  value  of  used  in  the  computation  for 

each  turbulence  level  as  computed  from  the  phase  structure  function.  The 

-2 

beam  wander  is  given  as  the  e  diameter  of  the  beam  position  distribution. 

Table  A-6  gives  the  transverse  coherence  length  and  the  re¬ 

sultant  long-term  beam  spread  0  .  Beam  spread  values  are  given  both  at 

-1  S  -2 
the  e  diameter  (0  )  (per  the  formula)  and  at  the  e  diameter  (0  ')  for 

s  s 

comparison  with  the  beam  wander.  We  have  computed  0g'  assuming  a  Gaussian 

distribution  (0  '  =  iJl  0  ).  In  order  to  make  the  result  more  general,  the 
s  s 

vacuum  beam  divergence  has  been  omitted. 

Table  A-7  summarizes  the  degradations  to  be  expected  in  the  illumi¬ 
nating  beam.  For  scintillation  we  include  the  point  intensity  fluctuation 
level  (dB(pt)),the  resolution  element  averaged  fluctuation  level  (dB(avg)), 
and  the  amplitude  correlation  length  p  .  Also  shown  are  the  short-term 
beam  wander  0,  the  short-term  beam  spread  y,  a^gy^he  long-term  beam 

spread  0  ' ,  each  given  at  the  e  ^  diameter  of  the  distribution  (y  is 
s 

the  root  square  difference  between  0  '  and  $) . 

For  scintillation  in  the  imaging  leg  (as  well  as  for  all  other  tur¬ 
bulence  effects  in  the  imaging  leg),  we  use  spherical  wave  formulae. 

2 

Table  A-8  shows  the  scintillation  results,  including  the  theoretical  a  , 

2 

the  value  corrected  for  saturation  a  (m) ,  the  value  corrected  for  source 
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Table  A-5.  Short-Term  Beam  Wander  over  66-km 
Slant  Path,  Free  Atmosphere 
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Table  A-6.  Long-Term  Beam  Spread  over  66-km 
Slant  Path,  Free  Atmosphere 
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averaging  a  (s) ,  and  the  receiver  averaged  value  a  (A)  (for  a  30-cm 
aperture).  Also  shown  are  the  amplitude  correlation  length  and  the 

resultant  rms  intensity  fluctuation  levels,  both  at  a  point  in  the  receiver 
plane  (dB(pt))  and  at  a  point  in  the  image  (dB(avg)). 

Table  A-9  gives  the  rms  image  dancing  computed  for  a  one-meter 

_2 

receiver  diameter,  both  at  the  rms  radius  (ct )  and  at  the  e  diameter  (a'). 
Note  that  this  parameter  is  independent  of  wavelength. 

To  calculate  long-term  point  spread  functions,  a  f  =  1.0  m,  D  = 

30  cm  diffraction-limited  imaging  lens  was  chosen.  Both  the  diffraction 

limit  and  the  turbulence  degraded  PSF's  were  computed.  Figure  A-12  shows 

typical  curves  (the  curves  are  normalized  to  unity  peak  intensity) .  From 

-2 

these  curves  and  the  accompanying  calculator  output,  the  e  diameter  of 

the  diffraction  limit  (0  )  and  the  degraded  (0.)  PSF's  were  obtained.  These 

O  A 

results  are  shown  in  Table  A-10,  along  with  the  atmospheric  degradation  0^ 

(root  square  difference  between  6^  and  0q)  and  the  on-axis  intensity 

reduction  I/I  .  To  evaluate  the  effect  of  lens  diameter  ( f //  constant)  on 
o 

image  degradation,  the  PSF  calculations  were  made  for  four  different 
diameters  (D  =  0.15,  0.3,  0.6,  and  0.9  m) .  Results  from  this  calculation 
are  shown  in  Table  A-ll  and  Figure  A-13. 

Table  A-12  summarizes  the  degradations  to  be  expected  in  the  image. 
For  scintillation  we  include  the  intensity  fluctuation  level  at  a  point  in 
the  receiver  plane  (dB(pt)),  the  fluctuation  level  at  an  image  point 
(dB(avg)),  and  the  amplitude  correlation  length  p  .  Also  shown  are  the 

X* 

image  dancing  a',  the  short-term  PSF  width  8.,  the  long-term  PSF  width 

A 

0^,  and  the  diffraction  limit  0q  (6  is  the  root  square  difference  between 

0  and  a').  Again,  these  results  are  for  a  f  =  1  m,  D  =  30  cm  imaging 

A  _2 
lens,  and  the  angular  quantities  are  given  at  the  e  diameter  of  the  re¬ 
spective  distributions. 
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Table  A-9.  Image  Dancing  over  66-km  Slant  Path, 
Free  Atmosphere 
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9.3 

Table  A-10.  Long-Term  Point  Spread  Function  Degradation 
over  66-km  Slant  Path,  Free  Atmosphere 
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Table  A-ll, 

Lens  Diameter  Dependence  of  Long- 

-Term  Point 

Spread 

Function 

Degradation 

over  66-km 

Slant 

Path,  Free  Atmosphere 

A  =  3.8  pm 

D 

80 

0A 

0D 

(cm) 

(prad) 

(prad) 

(prad) 

X/Io 
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10.9 
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60 
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0.564 

90  6.9  11.3  8.9  0.376 
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Figure  A-12.  Typical 
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Figure  A-13.  Lens  Diameter  Dependence  of  Point  Spread 
Function  Degradation  Over  66-km  Slant  Path, 

Free  Atmosphere 
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AIRCRAFT  BOUNDARY  LAYER  OPTICAL  TURBULENCE 

The  optical  degradations  introduced  by  propagation  through  the  air¬ 
craft  boundary  layer  will  be  qualitatively  similar  to  those  produced  by  the 
free  atmosphere.  We  are  primarily  interested  in  the  phase  degradations, 
i.e.,  the  wander  and  spread  in  the  illuminating  beam  and  the  dancing  and  reso¬ 
lution  loss  in  the  image.  Although  some  work  has  been  done  in  the  area  of 
optical  degradations  produced  bv  turbulent  boundary  layers  (Refs  34 ,  35), 
reliable  predictive  equations  are  not  generally  available.  We  have  there¬ 
fore  elected  to  utilize  experimental  results  to  obtain  upper  bound  estimates 
for  aircraft  boundary  layer  degradations. 

The  Air  Force  Weapons  Laboratory  (AFWL)  at  Kirtland  AFB  in  Albuquerque 

has  had  a  comprehensive  experimental  investigation  of  turbulence  degradations 

(including  both  boundary  layer  and  free  atmosphere  effects)  under  way  for 

some  time.  Although  the  major  findings  of  the  program  are  not  yet  published, 

useful  preliminary  results  were  made  available  to  the  SIS  program.  These 

2 

results  are  in  three  areas:  troposphere  C^  measurements  (discussed  pre¬ 
viously);  aircraft  boundary  layer  degradation  of  finite  beams;  and  boundary 
layer  MTF  measurements.  The  AFWL  beam  degradation  data  (only)  carry  a 
SECRET  security  classification. 

The  finite  beam  experiment  (Ref.  36)  involved  the  propagation  of 
a  10.6-micrometer  beam  from  the  Airborne  Laser  Laboratory  (ALL)  C-135  air¬ 
craft  to  a  target  T-39  aircraft  over  a  range  of  approximately  1  kilometer 
at  various  altitudes  and  air  speeds.  The  beam  was  focused  onto  a  target  board 
consisting  of  an  array  of  detectors  whose  temporal  output  was  recorded.  A 
simultaneous  recording  of  the  tracking  and  pointing  errors  in  the  ALL  C-135 
beam  projection  system  was  made.  Analysis  of  the  data  allowed  the  computation 
of  the  temporal  variation  of  both  the  beam  centroid  position  and  the  beam 
diameter  (as  well  as  other  parameters)  for  the  experimental  1-kilometer  path, 
and  subsequent  extrapolation  to  the  nominal  SIS  66-kilometer  slant  path. 

These  results  are  classified  SECRET,  to  be  consistent  with  the  AFWL  data 
upon  which  they  are  based.  To  avoid  classification  of  this  appendix,  those 
results  are  reporced  in  the  atmospheric  optical  effects  section  of  the  main 
body  of  this  repor*. 
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With  respect  to  imaging,  information  was  received  on  two  boundary 
layer  MTF  experiments  (Ref.  37),  both  designed  primarily  to  test  the  pre¬ 
dictions  of  simple  theory.  This  theory  predicts  that  the  MTF  has  the  form 


where  k  is  the  wave  number,  f  is  the  spatial  frequency  associated  with 

the  dominant  turbulence  scale,  and  a  is  a  strength  parameter  which  is 
proportional  to  the  dynamic  pressure.  The  line  spread  function  is  then 
predicted  to  show  saturation  at  apertures  large  compared  with  the  dominant 
scale  size  and  should  have  a  considerable  wavelength  dependence.  In  the 
first  experiment  a  0.6328  micrometer  beam  was  propagated  twice  through  the 
turbulent  boundary  layer  surrounding  a  specially  fitted  NASA  Lear  Jet  in 
flight.  Both  the  line  spread  function  (LSF)  and  the  MTF  were  measured  at 
a  variety  of  air  speeds  and  altitudes  and  for  beam  diameters  between 
10  and  24  mm.  In  some  measurements  an  aerodynamic  fence  was  placed  in  the 
alrstream  ahead  of  the  propagation  path.  The  LSF  degradation  as  measured  by 
the  peak  intensity  relative  to  the  non-degraded  peak  intensity  (l/lQ)  was 
between  0.73  and  0.86  micrometer  (one  way)  for  all  runs  without  the  fence. 
With  the  fence,  values  as  low  as  0.55  micrometer  were  reported.  The  ex¬ 
pected  dependence  on  dynamic  pressure  (altitude  and  air  speed)  was  present, 
only  when  the  fence  was  used.  The  LSF  saturation  with  aperture  size  was 
not  observed  even  though  the  24-mm  aperture  was  considerably  larger  than  the 
estimated  4.8-mm  dominant  scale  size.  The  fence,  often  used  to  quiesce 
flows  around  open  ports,  resulted  in  significantly  smaller  degradation. 

The  second  MTF  experiment  was  performed  in  a  wind  tunnel  and  was 
designed  to  investigate  further  the  wavelength  and  aperture  dependence. 
Although  analysis  is  not  complete,  preliminary  results  were  made  available 
which  verified  the  predicted  wavelength  scaling  (from  0.488  to  1.06  micro¬ 
meter)  and  which  showed  aperture  saturation  for  35  mm  beams. 
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In  order  to  use  the  results  from  these  experiments  in  obtaining 
estimates  of  boundary  layer  image  degradation,  we  assume  that  aperture  satura¬ 
tion  takes  place  at  all  wavelengths  before  40  mm.  (This  seems  reasonable 
since  boundary  layers  are  typically  5  cm  thick.)  Scaling  the  Lear  Jet  data 
to  40  mm,  degradations  of  0.45  should  be  expected  at  0.6328  micrometer. 

-  (a/A)2 

Using  a  1  -  e _ __  wavelength  dependence  (approximately  valid  for 

(a/A)2 

the  free  atmosphere),  we  compute  I/I  values  of  0.33,  0.73,  0.97,  and 
0.99  for  0.5145,  1.06,  3.8,  and  10.6  micrometers  respectively. 

COMBINED  RESULTS:  FREE  ATMOSPHERE  AND  AIRCRAFT  BOUNDARY  LAYER  TURBULENCE 

For  the  case  of  an  illumination  beam,  the  results  are  classified  and 
are  therefore  reported  in  the  body  of  this  report. 

For  the  imaging  case,  we  can  combine  the  boundary  layer  PSF  degrada¬ 
tion  by  taking  the  product  of  the  I/Iq  ratios  for  the  boundary  layer  and  the 
free  atmsophere  and  then  computing  resultant  PSF  widths.  Table  A-13  shows  the 
total  estimated  image  degradation  for  a  f  =  1  m,  D  =  30  cm  lens,  including 
the  scintillation  level  in  the  image  plane  (dB(avg)),  the  image  dancing  <*T, 
the  short-term  PSF  width  8^,  the  long-term  PSF  width  0^,  and  the  lens 

diffraction  limit  0  . 

o 
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Table  A-13.  Combined  Image  Degradation 


K 

o 

(um) 

dB(avg) 

(yrad) 

(yrad) 

(yrad) 

(yrad) 

0.5145 

0.10 

9.6 

21 

23 

2.82 

IQ'13 

1.06 

0.14 

9.6 

11 

15 

5.31 

3.8 

0.  30 

9.6 

21 

23 

20.8 

10.6 

0.39 

9.6 

58 

59 

58.1 

0.5145 

0.09 

9.3 

20 

22 

2.82 

-14 

1.06 

0.13 

9.3 

11 

14 

5.81 

10 

3.8 

0.25 

9.3 

21 

23 

20.8 

10.6 

0.32 

9.3 

58 

59 

58.1 

0.5145 

0.09 

9.3 

20 

22 

2.82 

in-i5 

1.06 

0.13 

9.3 

11 

14 

5.81 

10 

3.8 

0.25 

9.3 

21 

23 

20.8 

10.6 

0.31 

9.3 

58 

59 

58.1 
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Section  A- 5 

ATMOSPHERIC  AEROSOL  ATTENUATION  OVER  66-KM  SLANT  PATH 
(Summary  of  current  data  and  analysis) 


ATTENUATION  COEFFICIENTS* 

Figure  A-14  is  a  graph  of  aerosol  attenuation  coefficients  for  AFCRL 
atmospheric  aerosol  models,  at  sea  level,  for  a  normalized  "clear"  atmosphere 
(for  which  the  standard  visibility  range  is  23  kilometers). 

The  solid  curve  is  plotted  from  tabulated  values  taken  fromthe  AFCRL 
"Optical  Properties  of  the  Atmosphere"  report  (Ref.  1),  and  is  consistent 
with  Figure  22  of  that  report.  That  curve  is  also  reflected  in  the  aerosol 
coefficient  tabulations  of  other  reports  cf  the  1972  era  (e.g.,  AFCRL-72- 
0312  (Ref.  7)  and  AFCRL-72-061 1  (Ref.  4)  )  and  in  the  atmospheric  trans¬ 
mittance  computer  program  LOWTRAN  2,  AFCRL-72-0745  (Ref.  2). 

The  dotted  curve  marked  by  x's  was  plotted  from  tabled  values  (Ref.  19) 
supplied  to  SRL  by  AFCRL.  That  curve  represents  AFCRL' s  revised  model  for  an 
"average  continental  aerosol,"  based  upon  a  great  deal  of  experimental 
measurements  and  analysis  during  the  1972-74  time  period.  The  revised  data 
are  reflected  in  Figure  2  of  AFCRL-TR-74-0003  (Ref.  3)  but  are  not  yet  pub¬ 
lished  in  coefficient  tables  comparable  to  those  of  the  earlier  reports. 

They  do,  however,  appear  in  the  aerosol  spectral  data  listing  in  Computer 
Code  LOWTRAN  3  (Ref.  8),  recently  published  by  AFCRL. 

The  dotted  curve  of  Figure  A-14  marked  by  small  circles  represents 
an  AFCRL  estimate  of  an  "average  maritime  aerosol,"  according  to  tabulated 
values  (Ref.  20)  supplied  to  us  by  AFCRL.  The  data  are  the  "sea  level" 
values  applicable  to  the  "exchange  layer"  or  transition  region  between  the 
ocean  surface  "boundary  layer"  and  the  "upper  atmosphere,"  for  a  60%  salt 
spray,  40%  continental  aerosol  mixture.  Depending  upon  sea  state  and 
weather  (including  recent  weather  history),  the  exchange  layer  may  be  from 

*  Some  of  the  ideas  expressed  here  and  in  the  subsection  following  (entitled 
"Status  of  Models")  are  based  upon  personal  discussions  during  a  visit 
to  AFCRL  (see  Ref.  38). 
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1  to  5  kilometers  in  depth.  In  this  region,  the  atmospheric  aeronol  charac¬ 
teristics  are  considered  to  transition  from  a  predominantly  marine  type  in  the 
lower  part  to  a  predominantly  continental  type  at  the  top.  During  periods  of 
low  wind  and  a  relatively  calm  sea,  salt  particles  and  salt  water  droplets 
drift  downward  and  effectively  lower  the  top  of  the  marine  aerosol  region. 

The  plotted  curves  are  for  a  relative  humidity  of  75%.  However,  the 
variation  with  relative  humidity  is  negligible  for  continental  aerosols,  but 
is  significant  for  marine  aerosols.  This  is  due  to  the  hygroscopic  nature  of 
salt:  at  low  values  of  R.H.,  atmospheric  salt  content  consists  primarily  of 
solid  particles;  at  high  values  of  R.H. ,  much  of  the  atmospheric  salt  exists 
as  salt  water  droplets,  with  differing  optical  properties. 

STATUS  OF  MODELS 

The  plotted  (1974)  continental  aerosol  values  are  based  upon  many 
measurements  of  atmospheric  aerosol  particle  sizes,  chemical  determination  of 
constituent  materials,  and  measurement  of  the  complex  index  of  refraction  of 
those  materials.  Additionally,  optical  measurements  of  scattering  and 
attenuation  using  laser  sources  have  confirmed  the  accuracy  of  the  plotted 
results  at  representative  points  throughout  the  spectrum.  In  brief,  that 
curve  is  rather  well  established. 

in  contrast,  the  data  points  on  the  marine  aerosol  curve  are  cal¬ 
culated  values  based  upon  very  sketchy  data  about  marine  atmospheres.  Ex¬ 
perimental  verification  is  virtually  non-existent,  and  the  statistical 
variations  with  time,  altitude  and  location  are  almost  entirely  unknown.  In 
brief,  the  marine  aerosol  model  constitutes  a  "best  guess" — a  good  starting 
point  to  serve  as  a  basis  for  the  additional  research  and  experimentation 
required  to  establish  a  valid  data  base. 

APPLICATION  TO  SIS  SLANT  PATH  AEROSOL  PARTIAL  TRANSMITTANCE  ESTIMATES 

In  view  of  the  status  of  the  current  marine  aerosol  model,  the  reader 
is  cautioned  to  consider  the  following  results  as  rather  preliminary  (but 
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the  best  that  can  be  done  at  this  time).  The  cited  "average  continental" 
and  "estimated  marine"  models  were  employed  to  calculate  the  results  reported 
briefly  below  and  on  the  accompanying  graphs. 


Figure  A-15  is  a  graph  of  aerosol  transmittance  versus  wavelength 
for  a  66-kilometer  slant  path,  based  upon  the  revised  aerosol  coefficients 
of  Figure  A-14,  calculated  for  a  normalized  "clear”  atmosphere,  for  which 
the  sea  level  visibility  range  is  23  kilometers  (for  each  model).  The  solid 
curve  is  for  the  revised  "average  continental  aerosol  model"  and  the  dotted 
curve  assumes  that  the  estimated  marine  aerosol  model  (60%  sea  spray,  40% 
continental  aerosols)  exists  over  the  whole  path,  which,  from  our  conversations 
with  AFCRL, is  an  unlikely  situation.  The  dashed  curve  is  for  a  "composite" 
situation  in  which  the.  marine  aerosol  model  is  assumed  for  the  first  kilometer 
of  altitude  and  the  average  continental  model  is  assumed  for  the  remainder 
of  the  path.  This  is  considered  (postulated,  really)  to  represent  approxi¬ 
mately  a  "worst  case"  with  respect  to  atmospheric  aerosol  content. 

« 

Figure  A-16  is  a  similar  presentation  for  a  normalized  "hazy"  atmos¬ 
phere,  for  which  the  sea  level  visibility  range  is  5  kilometers  for  each 
model.  The  sea  level  aerosol  concentration  is  4.87  times  the  "clear"  value, 
with  a  decreasing  ratio  between  the  two  as  altitude  increases.  Above  5 
kilometers,  the  two  are  identical  in  the  models. 

Figure  A-17  is  a  similar  presentation  for  SRL's  normalized  "light 
haze"  atmosphere,  for  which  the  sea  level  visibility  is  10.8  km.  The  sea 
level  aerosol  concentration  for  this  visibility  is  the  geometric  mean  of  the 
"clear"  and  "hazy"  values,  or  2.2  times  the  "clear"  value  (2.2  =  x  1). 

Figures  A-18,  19,  and  20  present  the  same  curves  as  in  Figures  A-15, 

16,  and  17,  but  the  curves  are  grouped  to  compare  "clear",  "light  haze"  and 
"hazy"  transmittance  for  each  aerosol  type  in  lieu  of  the  previous  presentation 
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METHOD  OF  COMPUTATION 

Since  altitude  variation  tables  are  not  currently  available  for  the 
revised  continental  and  new  marine  models,  it  was  necessary  to  synthesize 
the  required  "optical  depth"  values  from  which  to  calculate  aerosol  trans¬ 
mittance  values  for  the  slant  path.  This  was  accomplished  in  the  following 
manner,  based  upon  the  essential  linearity  of  the  aerosol  attenuation 
process : 

(1)  It  was  determined,  from  data  tables  in  AFCRL-72-0497  (Ref.  1), 
that  the  equivalent  sea  level  path  length  (ESLPL)  of  a  0-1  km  vertical  path 
for  a  "clear"  atmosphere  is  0.664  km,  and  the  ESLPL  of  a  1-12  km  vertical 
path  is  0.656  km.  Similarly,  for  a  normalized  "hazy"  atmosphere  the  ESLPL 
for  0-1  km  is  0.604  km,  and  for  1-12  km  it  is  0.373  km. 

(2)  Then  for  our  66-kilometer  slant  path,  the  respective  ESLPL' s  are 
just  the  above  values  multipled  by  the  secant  of  the  path  zenith  angle 
(secant  79.35°  =  5.41).  Hence,  for  an  all-continental  or  all-marine  atmos¬ 
phere, *the  slant  path  optical  depth  values  are,  for  a  "clear"  atmosphere: 

OD  =  (0.664  +  0.656)5.41  y  =  7.14  y 

cl  Si 

For  a  normalized  "hazy"  atmosphere,  the  slant  path  "optical  depth"  is: 

OD  =  (0.604  +  0.373) (4.87) (5.41)  y  =  25./  y  , 

a  a 

where  y  is  the  "clear"  sea  level  value, 
a 

(3)  For  the  "composite"  atmosphere  the  two  assumed  layers  are  con¬ 
sidered  separately,  according  to  the  following: 

For  normalized  "clear": 

OD  =  (0.664  y  M  +  0.656  y  C)5.41 

3.  3 

I 

i 
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M  C 

where  y  and  y  are  the  sea  level  coefficient  values  for  the  marine  and 

cl  cl 

continental  aerosols,  respectively. 

For  normalized  "hazy": 

OD  =  [(0.604  x  4.87  y  M)  +  (0.373  x  4.87  y  C) ] 5 . 41 

a  a 

=  (2.94  y  M  +  1.817  y  C)5.41 
a  'a 

(4)  For  a  "light  haze"  atmosphere,  the  optical  depth  at  each  wave¬ 
length  was  taken  as  the  geometric  mean  of  the  "clear"  and  "hazy"  values.  This 
is  consistent  with  a  sea  level  aerosol  concentration  of  2.2  times  the  "clear" 
value,  where  the  hazy  value  is  4.87  times  the  clear  value,  with  comparable 
scaling  for  each  altitude  increment. 
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Figure  A-14.  Aerosol  Attenuation  Coefficients  (Average  Atmosphere  Models)  (Clear) 
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Aerosol  Partial  Transmittance  of  Path 


Figure  A-18.  Partial  Transmittance  of  Atmospheric  Aerosols, 
0  to  40,000  ft,  35  nmi  Slant  Path, 

Average  Continental  Aerosol  Model 
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Aerosol  Partial  Transmittance  of  Path 


Figure  A-20.  Partial  Transmittance  of  Atmospheric  Aerosols, 
0  to  40,000  ft,  35  nmi  Slant  Path, 

Composite  Aerosol  Model 
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Section  A- 6 


COMPUTER- CALCULATED  BAND  TRANSMITTANCE  CURVES 
FOR  SIS  66-KILOMETER  SLANT  PATH 

INTRODUCTION 

This  section  consists  of  70  computer-generated  plots  of  atmospheric 
transmittance,  calculated  for  the  SIS  nominal  66-kilometer  slant  path.  The 
plots  are  arranged  in  six  sets,  of  which  the  first  three  utilize  AFCRL's 
"average  continental  aerosol  model"  of  1974,  and  the  other  three  utilize 
AFCRL's  "estimated  marine  aerosol  model"  of  1975  (60%  sea  spray).  One  of 
the  three  wavelength  regions  is  covered  in  each  set,  with  the  wavelength 
regions  designated  as  Visible  and  Near  IR,  Mid-Range  IR,  and  Long-Wave  IR. 

A  wide  spectrum  of  atmospheric  conditions  is  represented  by  the  use  of 
AFCRL's  five  latitude-season  models  and  two  atmospheric  haze  models  ("clear) 
and  "hazy"),  with  an  intermediate  haze  level  ("light  haze")  also  included 
in  the  Visible  and  Near  IR  wavelength  region. 

The  "marine  aerosol"  sets  were  included  to  show  the  possible  effects 
of  sea  spray  in  the  atmosphere.  However,  the  uncertainties  of  the  marine 
aerosol  model,  as  discussed  previously,  dictate  cautious  use  of  the  marine 
results.  Also,  the  marine  aerosol  plots  assume  a  60%  sea  spray  aerosol 
throughout,  likely  a  very  pessimistic  assumption.  (The  relative  influence 
of  a  "composite"  aerosol,  with  sea  spray  only  in  a  lower  layer,  was  shown 
previously. ) 

These  LOWTRAN  plots  may  be  used  to  get  a  general  view  of  band  trans¬ 
mittance  over  the  wavelength  region  depicted,  and  the  variation  with  atmos¬ 
pheric  model  and  aerosol  (haze)  concentration.  Sets  I  and  IV  can  be  used, 
also,  to  estimate  laser  transmittance  at  visible  wavelengths  and  at  1.06 
micrometer  (because  of  the  absence  of  strong  molecular  absorption  lines  in 
the  vicinity).  However,  the  plots  will  not  give  accurate  estimates  for  DF 
laser  lines,  CO^  lines,  or  other  laser  frequencies  in  regions  having  strong 
molecular  resonance  absorption  lines.  The  reason  is  that  LOWTRAN  is  a  band 
transmittance  model;  it  computes  the  average  transmittance  over  an  interval 
of  20  wave-numbers.  The  transmittance  of  a  particular  laser  line  having 
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a  frequency  within  that  interval  may  be  somewhat  above  or  from  a  little  to 
far  below  this  average,  depending  upon  the  position  of  its  frequency  on  a 
high-resolution  plot  of  the  spectral  transmittance  in  the  interval.  Such 
monochromatic  transmittances  may  be  readily  calculated  manually  for  those 
laser  lines  for  which  accurate  attenuation  coefficients  have  been  computed 
and  validated. 

NOTE:  A  recent  research  report  (Ref.  39)  revises  the  previous 
estimates  of  water  vapor  continuum  absorption  in  the  8  to  13  micrometer 
wavelength  region,  and  includes  a  temperature  dependence  not  previously 
established.  The  new  data  was  incorporated  into  the  LOWTRAN  3  program, 
for  some  computer-calculated  runs  on  another  AFAL  program,  and  significant 
differences  in  the  results  were  found.  It  was  then  decided  to  rerun  the  pre¬ 
viously  generated  SIS  LOWTRAN  3  curves  (of  August-September  1975)  with  the 
cited  program  changes,  and  the  graphs  of  this  appendix  are  based  upon  the 
resulting  "Modified  LOWTRAN  3"  computer  program. 

The  curves  were  run  on  the  ASD  central  computer  facility  during 
April  1976,  using  the  terminal  in  Building  22B  and  the  Cal-Comp  plotter  in 
Building  22.  The  graphs  reproduced  herein  were  reduced  50%  in  size  from 
the  originals  to  facilitate  publication. 

The  wavelength  region  of  interest  was  divided  into  three  bands  for 
these  runs,  with  the  following  intervals  covered:  (1)  0.3  to  1.85  micrometers, 
(2)  1.85  to  5.5  micrometers,  and  (3)  7.1  to  14.3  micrometers.  The  wavelength 
region  from  5.5  to  /.I  micrometers  was  not  computed  because  it  is  known  to  be 
opaque  for  low  altitude  atmospheric  distances  of  more  than  a  few  hundred 
meters. 


The  above  cited  wavelength  intervals  are  approximate.  The  LOWTRAN 
computer  code  used  employs  frequency  (wave  numbers  per  centimeter)  rather 
than  wavelength,  and  computes  the  average  transmittance  over  20  wave-number 
intervals.  The  specific  frequency  regions  employed  for  the  three  cited 
bands,  and  the  computation  intervals  employed,  were  as  follows: 
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1.  5400  to  34,000  cm  (1.852  to  0.294  micrometers) 

Av  =  20  cm  \  5400  cm  ^  to  15,000  cm  ^ 

Av  =  200  cm  \  15,000  cm  ^  to  34,000  cm  ^ 

2.  1800  to  5400  cm  ^  (5.555  to  1.852  micrometers) 

Av  =  10  cm  ^ 

3.  700  to  1400  cm  ^  (14,286  to  7.143  micrometers) 

,  -1 
av  =  5  cm 

The  differing  computation  frequency  intervals  were  chosen  to  minimize 
computer  and  printout  time  while  preserving  adequate  resolution  in  the  graphs. 

For  each  wavelength  region  there  are  five  pages  of  graphs  for  each 
aerosol  type  (continental  and  marine).  Each  page  shows  the  transmittance 
plots  for  one  of  five  latitude-season  atmosphere  models,  at  each  of  two  or 
three  haze  levels,  keyed  to  visibility  range.  For  the  visible/near  IR 
wavelength  region  three  normalized  haze  levels  ("clear,"  "light  haze"  and 
"hazy")  are  shown,  for  which  the  sea  level  visibility  ranges  are  23,  10.8, 
and  5  kilometers,  respectively.  For  the  other  two  wavelength  regions,  only 
"clear"  and  "hazy"  graphs  are  shown  because  the  transmittance  variation  with 
haze  level  is  relatively  small  in  these  regions.  The  specific  aerosol  dis¬ 
tributions  making  up  the  normalized  "clear"  and  "hazy"  atmospheric  haze 
models  are  described  in  Ref.  1.  The  aerosol  particle  concentration  in  SRL's 
"light  haze"  model  is  the  geometric  mean  of  the  other  two  concentrations,  so 
that  the  sea  level,  aerosol  concentrations  are  in  the  ratio  of  1,  2.2,  and  4.87 
respectively,  for  "clear,"  "light  haze"  and  "hazy"  atmospheres. 

LATITUDE- SEASON  MODELS 

The  model  atmospheres  used  as  the  basis  for  LOWTRAN  computations 
are  given  in  Ref.  1.  Sea  level  values  of  major  variables  having  a  signifi¬ 
cant  influence  on  the  LOWTRAN  results  are  reported  in  Table  A-14. 
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Table  A-14.  Principal  Variables  in  Model  Atmosphere 


Subarctic 

Winter 

Model 

Mid-Lat . 

Winter 

Model 

Subarctic 

Summer 

Model 

Mid-Lat. 

Summer 

Model 

Tropical 

Model 

Pressure  (mb) 

1013 

1018 

1010 

1013 

1013 

Temperature 

-15.9°C 

-0. 8°C 

14°C 

21°C 

27°C 

Density 

(gm/m3) 

1372 

1301 

1220 

1191 

1167 

Water  vapor 
(gm/m3) 

1.2 

3.5 

9.1 

14 

19 

Ozone 

(pgm/m3) 

41 

60 

49 

60 

56 

CO  2 
(ppm) 

330 

330 

330 

330 

330 

LOUT RAN  2,  LOWTRAN  3,  AND  MODIFIED  LOWTRAN  3 

LOWTRAN  2  (Ref.  2)  is  an  atmospheric  transmittance  computer  model 
that  has  been  widely  used  during  the  past  three  or  four  years.  In  August 
and  September  1975  LOWTRAN  2  was  employed  for  slant  path  transmittance  cal¬ 
culations  for  the  SIS  nominal  66-kilometer  slant  path.  In  so  doing,  the 
following  revisions  were  made  in  the  basic  program  and  data  file: 

1.  Values  in  the  spectral  data  table  were  changed  per  an  up¬ 
date  received  from  AFCRL  on  3  March  1975. 

2.  The  aerosol  attenuation  function  was  deleted  and  tabulated  co¬ 
efficient  values  for  AFCRL's  Average  Continental  Aerosol  Model  (1974,  Ref. 
19)  and  Estimated  Marine  Aerosol  Model  (1975,  Rer.  20)  were  substituted, 
with  provision  for  selecting  either  model  at  the  start  of  a  run.  [The 
Average  Continental  Aerosol  Model  is  now  included  in  LOWTRAN  3  (Ref.  8).] 

3.  Provisions  were  made  for  employing  several  computation  intervals 
on  a  single  run,  for  different  wavelength  intervals. 
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A  linear  interpretation  routine  was  employed  with  the  aerosol  table 
and  other  appropriate  changes  were  made  to  accomplish  items  2  and  3  above. 


In  making  the  runs  for  this  report  an  additional  change  was  incorpor 
ated  in  the  program:  The  ozone  spectral  data  table  of  LOWTRAN  2  was  revised 
to  include  changes  made  in  the  recently  published  report  on  LOWTRAN  3  (Ref. 8). 

It  is  understood,  from  discussion  with  the  principal  author/researcher 
and  the  LOWTRAN  computer  programmer  at  AFCRL ,  that  incorporation  of  the 
above-cited  changes  into  LOWTRAN  2  effectively  convert  it  to  LOWTRAN  3  for 
most  applications,  and  we  have  labelled  the  graphs  of  this  report  accordingly. 
(The  published  LOWTRAN  3  also  incorporates  numerous  other  changes,  mostly  to 
make  the  program  more  flexible  and  useful  for  a  wider  range  of  applications.) 

"Modified  LOWTRAN  3"  is  SRL's  term  for  a  program  that  results  from  a 
change  in  the  8-13  micrometer  water  vapor  continuum  absorption  computation 
in  LOWTRAN.  This  change  is  made  to  incorporate  recent  research  findings 
reported  in  a  technical  paper  (Ref.  39)  presented  at  the  St.  Louis  IRIS 
meeting  on  3  February  1976.  The  required  program  changes  were  obtained  by 
phone  from  Mr.  Jim  Chetwynd,  AFCRL  LOWTRAN  programmer,  on  5  February  1976, 
with  first-order  verification  by  SRL.  They  remain  to  be  officially  verified, 
however,  in  a  formal  AFCRL  document. 

COMMENTS  ON  RESULTS 

The  first  page  of  Set  I  graphs  shows  the  rapid  drop-off  of  trans¬ 
mittance  at  the  short  wavelength  end  of  the  visible  region  due  to  ozone 
continuum  absorption  and  molecular  scattering.  It  also  shows  the  influence 
of  strong  molecular  resonance  absorption  in  a  number  of  near-IR  wavelength 
regions  beginning  at  the  top  end  of  the  visible  region  at  0.691  micrometer. 

In  comparing  the  three  curves  of  that  page,  the  very  strong  influence  of 
atmospheric  haze  level  upon  visible  wavelengths  is  apparent,  as  is  the  de- 
creasing  effect  of  this  factor  in  the  near  XR* 

The  remaining  pages  of  Set  I  exhibit  only  small  changes,  showing 
the  relatively  negligible  influence  of  latitude  and  season  upon  visible  and 
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near  IR  wavelengths.  (Note  the  difference  in  amplitude  of  some  of  the 
structure,  however;  e.g.,  the.  spike  at  1.45  micrometer,  and  the  width  of 
the  1.4  to  1.8  micrometer  window.) 

The  first  page  of  Set  II  shows,  for  the  66-kilometer  path,  the  three 
prominent  windows  in  the  mid-IR  region.  Comparison  of  the  three  curves  shows 
the  reduced  influence  of  atmospheric  haze  in  this  wavelength  region,  as  com¬ 
pared  to  the  visible  and  near  IR  region  of  the  previous  graphs. 

The  remaining  four  pages  of  Set  II  show  substantial  changes,  illus¬ 
trating  the  relatively  greater  influence  of  latitude-season  variables  upon 
atmospheric  transmittance  in  this  region.  The  amount  of  atmospheric  water 
vapor  (absolute  humidity)  accounts  for  the  major  differences  in  results  be¬ 
tween  the  five  latitude-season  models,  with  the  largest  difference  being  at 
the  longer  wavelengths.  The  relatively  good  transmittance  of  the  2.1-2. 3 
micrometer  and  3.5-4  micrometer  atmospheric  windows,  over  a  wide  range  of 
atmospheric  variables,  is  an  important  feature  shown  by  the  five  pages  of 
Set  II. 


Set  III  shows  the  wide  IR  window  from  about  8  to  about  14  micrometers, 
most  often  used  for  passive  thermal  imaging.  The  ozone  absorption  band  be¬ 
tween  9  and  10  micrometers  is  evident;  however,  the  major  absorber  is  water 
vapor,  with  CO^  also  significant.  The  two  curves  of  each  page  show  the 
relatively  small  influence  of  atmospheric  haze  upon  this  wavelength  region. 

The  five  pages  of  Set  III  illustrate  the  large  influence  of  latitude- 
season  variables  in  this  region.  The  chief  constituent  responsible  for  the 
large  variation  is  water  vapor  content.  Because  of  this  influence,  and  the 
large  dependence  of  atmospheric  water  vapor  capacity  upon  air  temperature, 
high  values  of  transmittance  are  most  likely  in  cold  weather  and  low  values 
are  most  likely  in  warm  weather.  The  temperature  itself  has  an  influence 
upon  the  optical  absorption  of  molecular  constituents,  as  evidenced  in 
several  references  (e.g..  Refs  14  and  20),  and  this  influence  is  reflected 
in  the  revised  water  vapor  continuum  absorption  incorporated  into  "Modified 
LOWTRAN  3." 
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Sets  IV,  V,  and  VI  duplicate  Sets  I,  II,  and  III,  respectively, 
but  with  the  AFCRL  Estimated  Marine  Aerosol  Model  (Ref.  20)  data  table 
substituted  for  the  continental  aerosol  model  table.  The  major  influence 
of  the  sea  spray  in  the  marine  model  is  seen  to  be  a  much  reduced  average 
transmittance  in  the  near  and  mid-IR  regions,  as  one  would  predict  from  the 
data  of  Section  V  of  this  report.  The  reader  is  again  cautioned,  however, 
of  the  uncertainties  in  the  marine  aerosol  model  and  the  likely  great  varia¬ 
bility  in  over-ocean  aerosols  (as  compared  to  continental  aerosols.)  resulting 
from  sea  state  and  weather  variables.  The  marine  aerosol  LOWTRAN  curves 
may  present  an  overly-pessimistic  picture  except  for  the  most  severe  conditions. 
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Set  I  (cont.) 


i 

Visible  and  Near  IR  Computer-Calculated  Band  Transmittance  over  66-km  Slant 
Path,  Mid-Latitude  Summer  Model,  Average  Continental  Aerosol  Model. 
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Visible  and  Near  IR  Computer-Calculated  Band  Transmittance  over  66-km  Slart 
Path,  Sub-Arctic  Summer  Model,  Average  Continental  Aerosol  Model. 
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Set  I  (cont.) 


Visible  and  Near  IR  Computer-Calculated  Band  Transmittance  over  66-km  Slant 
Path,  Mid-Latitude  Winter  Model.  Average  Continental  Aerosol  Model. 
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MRVElENOtM 


Visible  and  Near  IP  Computer-Calculated  Band  Transmittance  over  66-km  Slant 
Path,  Sub-Arctic  Winter  Model,  Average  Continental  Aerosol  Model. 


Mid-Range  IR  Computer-Calculated  band  Transmittance  over  66-km  Slant  Path, 
Mid-Latitude  Summer  Model,  Average  Continental  Aerosol  Model. 
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Mid-Range  IR  Computer-Calculated  Band  Transmittance  over  66-km  Slant 
Sub-Arctic  Winter  Model,  Average  Continental  Aerosol  Model. 
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Long-Wave  IR  Computer-Calculated  Band  Transmittance  over  66-km  Slant  Path, 
Tropical  Atmosphere  Model,  Average  Continental  Aerosol  Model. 
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Set  III  (cont.) 


Long-Wave  IR  Computer-Calculated  Band  Transmittance  over  66-km  Slant  Path 
Mid-LatituSe  Summer  Atmosphere  Model,  Average  Continental  Aerosol  Model. 


m 


Long-Wave  IR  Computer-Calculated  Band  Transmittance  over  66-km  Slant  Path, 
Sub-Arctic  Winter  Model,  Average  Continental  Aerosol  Model. 
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Set  IV  (LOW IRAN  3  -  Marine) 


vaIIABIF 


WRVEIENGTH 

LIGHT  HAZE 

_  m  1 

WRVEl ENO’H 


^  10  0  36  0.60  0.66  0  >0  0  >6  l.!0  1.26  1  <0  I  66  1  70  1  66  2.00 

KRVEIENOTH 

Visible  and  Near-IR  Computer-Calculated  Band  Transmittance  over  66-km  Slant 
Path,  Tropical  Atmosphere  Model,  Estimated  Marine  Aerosol  Model. 
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Set  IV  (cont.) 


MRVEIENOTH 


Visible  and  Near-IR  Computer-Calculated  Rand  Transmittance  over  66-km  Slant 
Path,  Mid-Latitude  Summer  Model,  Estimated  Marine  Aerosol  Model. 
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Set  IV  (cont.) 


Visible  and  Near-IR  Computer-Calculated  Band  Transmittance  over  66-km  Slant 
Path.  Sub-Arctic  "Summer  Model,  Estimated  Marine  Aerosol  Model. 


Set  IV  (cont.) 
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Con,pucer-CalcUUted  Band  Transmittance  over  66-km  Slant 
Pdth,  Mid-Latitude  Winter  Model,  Estimated  Marine  A.rn.m  „odel,  *“ 
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Set  IV  (cont.) 
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WRVELENOTH 


i  ns  7  QO 


Visible  and  Near-IR  Computer-Calculated  Band  Transmittance  over  66-km  Slant 
Path,  Sub-Arctic  Winter  Model,  Estimated  Marine  Aerosol  Model. 
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Mid-Range  IR  Computer-Calculated  Band  Transmittance  over  66-km  Slant  Path 
Mid-Latitude  Summer  Model,  Estimated  Marine  Aerosol  Model. 
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Set  V  (cent.) 


Mid-Range  IR  Computer-Calculated  Band  Transmittance  over  66-ktn  Slant  Path, 
Sub-Arctic  Summer  Model,  Estimated  Marine  Aerosol  Model. 
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1.70 


HRVElENOrH 


^11d^Ran£e_I£  Computer-Calculated  Band  Transmittance  over  66-km  Slant  Path 
Sub-Arctic  Winter  Model,  Estimated  Marine  Aerosol  Model. 
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Set  VI  (Modified  LOWTRAN  3  -  Marine) 
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wnvfifNOtH 


URVELENOTH 


Lon^Wave^IR  Computer-Calculated  Band  Transmittance  over  66-km  Slant  Path, 
Tropical  Atmosphere  Model,  Estimated  Marine  Aerosol  Model. 
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Long- Wave  IR  Computer-Calculated  Band  Transmittance  over  66-km  Slant  Path, 
Sub-Arctic  Summer  Model,  Estimated  Marine  Aerosol  Model. 
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BtSrWWV  ABLE  COPY 

Set  VI  (cont . ) 

J 


LongMteve^Ijl  Computer-Calculated  Band  Transmittance  over  66-km  Slant  Path, 
Mid-Latitude  Winter  Model,  Estimated  Marine  Aerosol  Model. 
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Section  A- 7 


COMPARISON  OF  MEASURED  AND  CALCULATED  DF  LASER  LINE  TRANSMITTANCES 

INTRODUCTION 

Comparisons  were  made  of  a  set  of  DF  laser  line  atmospheric  trans¬ 
mittance  measurements,  with  several  sets  of  calculated  values.  The  atmos¬ 
pheric  water  vapor  content  during  the  measurements  was  nearly  equal  to  that 
of  the  AFCRL  mid-latitude  summer  model  atmosphere;  hence  comparison  with  that 
model  should  be  valid  and  significant.  The  measurements  were  made  by  NRL 
over  a  5-kilometer  sea  level  path  (over  water)  at  Cape  Canaveral  in  March 
of  1975,  and  the  results  (Ref.  40)  were  furnished  SRL  in  the  form  of  a 
graph  of  attenuation  coefficients  (Figure  A-21) ,  derived  from  the  experi¬ 
mental  transmittance  measurements.  Figure  A-21  also  shows  NRL-calculated 
molecular  absorption  coefficients,  and  these  were  supplied  in  tabular  form 
as  well. 


The  OSU  (Ref.  16)  and  AFCRL  (Ref.  41)  coefficients  used  were  from 
unpublished  data  supplied  to  us  on  8  August  1975  and  11  September  1975, 
respectively,  and  consist  of  computer-calculated  values  of  molecular  absorp¬ 
tion  coefficients  for  OSU-revised  line  frequencies.  These  data  encompass 
36*  DF  laser  lines,  which  include  the  22  lines  in  the  NRL  measurements. 

There  should  be  a  reasonably  constant  difference  between  the  experi¬ 
mental  and  calculated  curves  of  Figure  A-21,  attributable  to  aerosol  absorp¬ 
tion  and  scattering.  The  data  of  Figure  A-21  (including  the  "extinction" 
values  for  1.06  micrometer)  were  used  in  conjunction  with  "continental"  and 
"marine"  aerosol  models  (previously  described)  to  estimate  aerosol  attenua¬ 
tion  coefficient  values  for  the  measurements  program  atmosphere.  (It  would 
have  been  helpful  if  more  data  were  available  for  this  purpose;  e.g., 
"visibility  range"  or  extinction  measurements  at  a  frequency  in  the  visible 
spectrum  region,  and  time  correlation  of  the  DF  measurements  with  the  1.06 
micrometer  measurements.) 


*The  AFCRL  table  for  the  DF  laser  1-0  P(9)  is  erroneous;  hence  valic  data 
cover  35  lines. 
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Wavelength  (vim) 

Figure  A-21.  Preliminary  DF  Extinction  Data,  CCAFS ,  5.08  km 


RESULTS 


Perhaps  the  most  meaningful  comparison  of  measured  and  calculated 
values  is  that  of  Figure  A-22.  This  bar  graph  shows  the  following  values, 
respectively,  for  each  of  the  22  DF  laser  lines: 

1.  The  NRL-measured  transmittance  for  their  5-kilometer  sea  level 
path,  as  derived  from  coefficient  values  scaled  from  Figure  A-21. 

2.  A  transmittance  value  based  upon  NRL-calculated  molecular  absorp¬ 
tion  coefficients,  together  with  the  estimate  of  the  aerosol  attenuation 
coefficient  (y  )  for  the  test  path. 

cl 

3.  A  transmittance  value  based  upon  AFCRL-calculated  molecular 
absorption  coefficients  for  the  revised  OSU  line  frequencies,  together 
with  the  SRL  estimate  of  aerosol  attenuation  coefficient. 

The  results  are  shown  in  a  different  form  in  Figures  A-23  and  A-24, 
in  which  differences  between  measured  values  of  total  attenuation  coefficient 
(y)  and  calculated  molecular  absorption  coefficient  are  plotted.  Figure 
A-23  was  copied  from  an  NRL-supplied  graph.) 

COMMENTS  ON  RESULTS 

Referring  again  to  Figure  A-22,  there  appears  to  be  good  agreement 
(T  values  within  0.05)  between  measured  and  calculated  transmittances  for 
17  of  the  22  lines  when  the  AFCRL/OSU  coefficient  values  are  used.  [All 
except  1-0  P (6) ,  1-0  P (7)  ,  2-1  P(4),  1-0  P(8)  and  1-0  P (10)  .  ] 

Use  of  the  AFCRL/OSU  coefficients  in  lieu  of  the  NRL  values  results 
in  much  closer  agreement  between  calculated  and  measured  transmittances  for 
the  2-1  P[7)  and  2-1  P(10)  lines,  with  lesser  improvement  for  three  other 
lines,  and  small  degradation  for  three  lines. 

Assuming  that  the  measured  values  are  accurate,  that  our  estimate  of 
y  is  reasonable,  and  that  there  was  no  significant  unmeasured  and  unknown 
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Figure  A-22.  Measured  and  Calculated  Transmittances  of  22  DF  Laser  Lines 

Over  5-Kilometer  Sea  Level  Path 
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disparity  between  the  test  atmosphere  and  the  mid-latitude  summer  model, 
the  following  conclusions  are  postulated  with  regard  to  calculations  of  DF 
laser  atmospheric  transmittance  using  the  AFCRL/OSU  coefficients. 

1.  Calculated  transmittances  for  1-0  P ( 8 )  and  1-0  P(10)  may  be 
significantly  high. 

2.  Calculated  transmittances  for  1-0  P(6),  1-0  P(7),  and  2-1  P(4) 
may  be  significantly  low. 

In  each  case  where  there  is  a  significant  disparity  between  measured 
and  calculated  values,  the  laser  line  frequency  lies  on  or  near  a  steep 
slope  of  a  molecular  absorption  line  in  the  atmospheric  transmission  spec¬ 
trum.  Hence  a  very  small  error  in  the  laser  line  frequency  or  in  the  spec¬ 
tral  characteristics  of  the  absorption  line  can  make  a  significant  difference 
in  the  calculated  transmittance  in  these  regions.  Figures  A-25  through  A-29 
are  high  resolution  transmittance  spectra  (from  Ref.  3)  which  illustrate 
this  phenomenon.  Each  of  the  36  DF  laser  lines  for  which  OSU-revised  line, 
frequencies  are  available  is  marked  on  these  graphs. 

On  the  foregoing  basis,  SRL  concludes  that  the  recently  computed 
OSU  or  AFCRL  absorption  coefficients  (based  on  the  OSU-measured  line 
frequencies)  provide  an  adequate  (and  probably  best  currently  available) 
basis  for  estimating  the  atmospheric  transmittance  (molecular  absorption 
component)  of  DF  lasers.  Estimates  will  likely  be  more  accurate  for  laser 
lines  whose  frequencies  are  significantly  distant  (e.g.,  a  tenth  or  so  of 
a  wave  number)  from  steep  slopes  in  the  atmospheric  transmittance  spectrum. 

It  would  be  well  to  make  this  factor  a  consideration  in  selection  of  laser 
lines  for  an  operational  DF  laser,  in  addition  to  the  estimated  atmospheric 
transmittance  (over  a  range  of  conditions)  and  other  factors  involved  in 
overall  system  performance  and  calculation  thereof. 

METHOD  OF  CALCULATION 

Measured  attenuation  coefficient  values  (y)  were  scaled  directly 
from  Figure  A-21,  and  converted  to  "measured"  5-kilometer  path  transmittances 
by  the  expression 


calculated  value  using  revised  f's  and  coefficients 
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Figure  A-25.  Atmospheric  Transmittance  Due  to  Molecular  Absorption 
Through  a  10-kra  Horizontal  Path  at  Sea  Level 


3P10 


Figure  A-26.  Atmospheric  Transmittance 
Through  a  10-km  Horizontal 


ulated  value,  using  revised  f  and 


Figure  A  28.  Atmospheric  Transmittance  Due  to  Molecular  Absorption  Through 

a  10-km  Horizontal  Path  at  Sea  Level 


OSU-calculated  and  NRL-calculated  values  of  molecular  absorption 
coefficients  were  taken  directly  from  tabled  values  supplied  by  OSU  (Ref.  16) 
and  NRL  (Ref  40).  Since  the  AFCRL  values  did  not  include  water  vapor  con¬ 
tinuum,  OSU-reported  values  of  this  coefficient  were  added  to  the  AFCRL 
molecular  absorption  coefficients. 

Tables  A-15  and  A-16  summarize  the  calculations  supporting  the 
graphical  presentation  of  results.  Table  A-17  is  a  copy  of  the  OSU  laser 
line  frequencies  and  coefficients  compilation.  Table  A-18  is  a  tabulation 
of  molecular  absorption  coefficients  by  altitude  layers,  for  the  2-1  P(7) 

DF  laser  line.  This  is  one  of  36  such  tables  furnished  by  AFCRL  (Ref.  41). 
These  data  replace  the  corresponding  molecular  absorption  coefficient  data 
of  Reference  7. 


Table  A-15.  Measured  Attenuation  Coefficients  and  Three 
Calculated  Molecular  Absorption  Coefficients  of  22  DF 
Laser  Lines  over  5-km  SL  Path 


y*  3*C 

Y  Y*  Y*  AFCRL* 


Line 

Id 

NRL 

(meas' d) 

NRL 

(calc' d) 

Ay 

NRL 

OSU 

(calc ' d) 

Ay 

OSU 

Revised 
(calc' d) 

ay 

AFCRL 

3-2  P(10) 

0.076 

0.0532 

0.0228 

0.0491 

0.0269 

0.0503 

0.0257 

3-2  P (9) 

0.069 

0.0354 

0.0336 

0.0357 

0.0333 

0.0359 

0.0331 

2-1  P  (12) 

0.063 

0.0337 

0.0293 

0.0351 

0.0279 

0.0340 

0.0290 

3-2  P (8) 

0.091 

0.0514 

0.0396 

0.0504 

0.0406 

0.0450 

0.0460 

2-1  P(ll) 

0.066 

0.0351 

0.0309 

0.0365 

0.0295 

0.0323 

0.0337 

3-2  P  (7) 

0.101 

0.0604 

0.0406 

0.0646 

0.0364 

0.0665 

0.0345 

2-1  P(10) 

0.110 

0.0468 

0.0632 

0.0693 

0.0407 

0.0727 

0.0373 

3-2  P (6)  , 

0.066 

0.0288 

0.0372 

0.0292 

0.0368 

0.0289 

0.0371 

2-1  P (9) 

0.083 

0.0415 

0.0415 

0.0397 

0.0433 

0.0383 

0.0447 

3-2  P (5) 

0.050 

0.0222 

0.0278 

0.0216 

0.0284 

0.0216 

0.0284 

2-1  P (8) 

0.056 

0.0327 

0.0233 

0.0344 

0.0216 

0.0345 

0.0215 

2-1  P (7) 

0.134 

0.0645 

0.0695 

0.0935 

0.0405 

0.0907 

0.0433 

1-0  P(10) 

0.094 

0.0366 

0.0574 

0.0371 

0.0569 

0.0370 

0.0570 

2-1  P (6) 

0.099 

0.0707 

0.0283 

0.0596 

0.0394 

0.0580 

0.0410 

1-0  P(9) 

0.072 

0.0403 

0.0317 

0.0377 

0.0343 

2-1  P (5) 

0.061 

0.0282 

0.0328 

0.0286 

0.0324 

0.0288 

0.0322 

1-0  P(8) 

0.178 

0.1074 

0.0706 

0.115 

0.0630 

0.1102 

0.0678 

2-1  P(4) 

0.071 

0.0632 

0.0078 

0.0552 

0.0158 

0.0546 

0.0164 

1-0  P (7) 

0.054 

0.0533 

0.0007 

0.0604 

-0.0064 

0.0603 

-0.0063 

2-1  P ( 3) 

0.071 

0.0409 

0.0301 

0.0407 

0.0303 

0.0402 

0.0308 

1-0  P (6) 

0.096 

0.0849 

0.0111 

0.0874 

0.0086 

0.0833 

0.0127 

1-0  P (5) 

0.086 

0.0679 

0.0181 

0.0716 

0.0144 

0.0691 

0.0169 

*:  All  mid-latitude  summer,  sea  level. 
**:  OSU  values  for  1^0  continuum  added. 
Ay's  are  NRL-measured  (-)  calculated. 
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Table  A-16.  Measured  Transmittances  and  Two  Calculated 
Transmlttances  of  22  DF  Laser  Lines  over  5-km 
SL  Path  (Mid-latitude  Summer) 


T  T  T 


Line 

Id 

NRL** 

(meas ' d) 

Est '  d 

Ya 

NRL 

(calc'd) 

SRL* 

(calc'd) 

3-2  P(10) 

0.684 

0.034 

0.647 

0.656 

3-2  P (9) 

0.708 

0.034 

0.707 

0.  705 

3-1  P (12) 

0.730 

0.034 

0.713 

0.712 

3-2  P(8) 

0.634 

0.034 

0.652 

0.674 

2-1  P (11) 

0.719 

0.034 

0.708 

0.718 

3-2  P ( 7) 

0.604 

0.034 

0.624 

0.605 

2-1  P(10) 

0.577 

0.034 

0.668 

0.587 

3-2  P (6) 

0.719 

0.034 

0.731 

0.730 

2-1  P (9) 

0.660 

0.034 

0.686 

0.697 

3-2  P (5) 

0.779 

0.034 

0.755 

0.757 

2-1  P (8) 

0.756 

0.034 

0.716 

0.710 

2-1  P(7) 

0.512 

0.034 

0.611 

0.536 

1-0  P(10) 

0.625 

0.034 

0.703 

0.701 

2-1  P(6) 

0.610 

0.034 

0.592 

0.631 

1-0  P(9) 

0.700 

0.034 

0.690 

0.699 

2-1  P (5) 

0.737 

0.034 

0.733 

0.731 

1-0  P(8) 

0.411 

0.033 

0.496 

0.489 

2-1  P(4) 

0.  701 

0.033 

0.618 

0.645 

1-0  P(7) 

0.763 

0.032 

0.653 

0.630 

2-1  P (3) 

0.701 

0.031 

0.698 

0.700 

1-0  P(6) 

0.619 

0.030 

0.563 

0.556 

1-0  P (5) 

0.651 

0.029 

0.616 

0.612 

*:  Using  revised  OSU  line  frequencies  and  revised  AFCRL  coefficients 
(11  September  1975). 

**:  Based  upon  NRL  graph  of  measured  y's* 
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Table  A-17.  Mid-Latitude  Summer  Sea  Level  Absorption  Coeffients 


Table  A-18.  Molecular  Absorption  Coefficients  by  Altitude  Layers 

for  DF  Laser  2-1  P(7)  Line 
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Section  A-8 


COMPARISON  OF  MEASURED  AND  CALCULATED  CO  LASER 
LINE  TRANSMTTTANCES  L 


INTRODUCTION 

Comparisons  were  made  of  a  set  of  CO^  laser  line  atmospheric  trans¬ 
mittance  measurements  (Ref.  42)  5  with  several  sets  of  calculated  values. 

The  measurements  were  made  by  NRL  over  a  5-kilometer  sea  level  path  (over 
water)  at  Cape  Canaveral,  and  the  results  of  the  measurements  were  supplied 
to  SRL  in  the  form  of  a  graph  of  attenuation  coefficients  derived  from  the 
experimental  transmittance  values.  Calculated  values  for  part  of  the  lines 
were  also  included  on  the  graph. 

It  was  desired  to  compare  the  NRL-supplied  data  with  estimated  trans- 
mittances  for  the  same  path  calculated  for  AFCRL  model  atmospheres.  The 
NRL  graph  indicates  measured  path  water  vapor  partial  pressure  and  CO^ 
concentration  of  15.8  torr  and  350  ppm,  respectively,  while  NRL  calculations 

were  accomplished  for  17  torr  and  350  ppm.  The  AFCRL  mid-latitude  summer 

3  3 

model  has  14  gm/m  H^O  water  vapor  and  the  tropical  model  has  19  gm/m 

HO  vapor.  The  CO  concentration  is  fixed  at  330  ppm  in  the  AFCRL  models. 

^  ^  3 

A  partial  pressure  of  one  torr  is  approximately  equivalent  to  1  gm/m  H2O 

at  20°  C,  and  departs  but  little  from  that  value  over  a  fairly  wide  tempera¬ 
ture  range.  Then,  since  these  two  constituents  are  the  principal  absorbers 
in  this  wavelength  region,  the  transmittance  of  the  test  path  could  be  ex¬ 
pected  to  be  about  midway  between  that  of  the  mid-latitude  summer  and  tropical 
model  atmospheres. 

Model  atmosphere  molecular  absorption  coefficient  data  (Ref.  4)  were 
available  for  12  of  the  CO^  lines  for  which  NRL  measurements  were  made; 
hence  comparisons  were  possible  at  those  corresponding  wavelengths.  A 
nominal  value  of  aerosol  attenuation  coefficient  (0.013  per  km)  for  a 
normalized  clear  atmosphere  was  added  to  each  tabled  value  of  molecular 
absorption  coefficient,  based  upon  AFCRL' s  "average  continental  aerosol 
model"  (1974),  to  arrive  at  a  total  attenuation  coefficient.  Transmittances 
for  the  5-km  path  were  calculated  from  the  resulting  values  of  attenuation 
coefficient . 
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RESULTS 


The  results  are  graphed  on  Figure  A- 30.  The  first  bar  for  each  wave¬ 
length  shows  the  NRL-measured  transmittance,  and  the  second  bar  (where 
given)  shows  the  corresponding  NRL-calculated  transmittance,  both  derived 
from  the  NRL  graph  of  attenuation  coefficients.  Then,  where  given,  the  two 
horizontal  marks  adjacent  to  the  second  bar  give  the  SRL-calculated  trans- 
mi*"tances  for  the  AFCRL  mid-latitude  summer  model  (upper  mark)  and  AFCRL 
tropical  model  (lower  mark). 

Table  A-19  summarizes  the  calculations  and  tabulated  results  support¬ 
ing  the  graphical  presentation. 

COMMENTS  ON  RESULTS 

Computed  estimates  of  transmittance  for  the  P-40  line  are  greatly  in 
error  in  comparison  to  the  measured  value.  A  possible  explanation  may  lie 
in  the  close  proximity  of  a  strong  water  vapor  absorption  line  to  the  P40 
frequency.  (See  Figure  A-31,  adapted  from  Ref.  4.)  A  small  error  in  the 
estimated  frequency  of  either  line  could  result  in  near  superposition  of 
the  two  attenuation  spikes  shown  at  wave  numbers  of  925  and  925.4  per  centi¬ 
meter  on  Figure  A-31,  and  could  account  for  a  much  greater  attenuation  of 
the  P40  line  than  is  accounted  for  by  the  available  (AFCRL-computed) 
attenuation  coefficients. 

It  also  appears,  from  Figure  A- 30,  that  the  SRL-calculated  "bounds" 
(based  upon  the  AFCRL  tropical  and  mid-latitude  summer  models)  are  some¬ 
what  optimistic  for  the  P28  line,  somewhat  pessimistic  for  the  P10,  P14, 
and  P18  lines,  and  "about  right"  for  the  remaining  lines  for  which  atten¬ 
uation  coefficients  were  available. 

The  following  observations  are  made  comparing  these  results  and 
comparable  results  for  a  group  of  DF  laser  wavelengths  in  the  neighborhood 
of  3,8  micrometer,  reported  in  Section  A-7  of  this  appendix. 
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1.  The  "limiting"  transmittance  (i.e.,  the  envelope  of  peak  trans¬ 
mittance  values  over  the  wavelength  region  of  interest)  is  much  lower  in 

the  CO^  laser  wavelength  region  than  it  is  in  the  DF  laser  wavelength  region, 
mainly  because  of  the  much  higher  water  vapor  continuum  absorption  in  the 
CO^  wavelength  region. 

2.  CC>2  laser  wavelengths  are  at  CC^  absorption  line  frequencies 
(see  Figure  A-31);  hence  the  transmittance  of  most  CO^  laser  lines  is  poorer 
than  the  average  "band  transmittance"  for  the  adjoining  wavelength  region. 

3.  The  wavelengths  of  DF  laser  lines  having  "good"  or  better  atmos¬ 
pheric  transmittance  do  not  coincide  with  atmospheric  absorption  lines;  hence 
the  transmittance  of  such  lines  is  better  than  the  average  "band  transmittance" 
for  the  adjoining  wavelength  region. 

4.  The  following  brief  summary  of  NRL  measured  transmittances  for 
the  5-kilometer  overwater  path  at  Cape  Canaveral  shows  one  comparison  of 
DF  laser  and  CO^  laser  atmospheric  transmittances,  for  roughly  comparable 
atmospheric  conditions: 

a.  For  all  22  DF  lines  employed,  measured  transmittances  fell  be¬ 
tween  0.41  and  0.79.  Twenty  were  above  0.6 

b.  For  all  16  CO2  lines  employed,  measured  transmittances  fell 
between  0.0002  and  0.25.  The  10.591  micrometer  value  was  0.11  and 
11  of  13  lines  between  10.2  micrometer  and  10.9  micrometer  were 
between  0.095  and  0.15. 
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Figure  A-30.  Measured  and  Calculated  Values  of  Transmittance 
Over  a  5-km  Over-Water  Path  for  a  Group  of  C0o  Laser  Lines 

12 


Figure  A-31.  Atmospheric  Transmittance  Due  to  Molecular  Absorption  Through 

a  10-km  Horizontal  Path  at  Sea  Level 


Table  A-19.  Measured  and  Calculated  Values  of  Attenuation 
Coefficient  and  Transmittance  over  5-km 
Over-Water  Path,  for  a  Group 
of  CO^  Laser  Lines 


Attenuation 

Coefficient,  y, 

and  Transmittance 

T _ 

CO  2 

SRL  Calc ' d 

SRL  Calc 

'd 

Laser 

X 

NRL  Measured 

NRL  Calc ' d 

(Mid-lat 

summer) 

(Tropical) 

Line 

(pm) 

v* 

T 

y* 

T 

Y** 

T 

y* 

T 

R30 

9.2197 

0.275 

0.253 

R20 

9.2714 

0.315 

0.207 

P20 

9.5524 

0.345 

0.178 

R18 

10. 2604 

0.395 

0.139 

0.410 

0.129 

0.338 

0.185 

0.535 

0.069 

P6 

10.4583 

0.380 

0.150 

0.358 

0.167 

0.323 

0.199 

0.526 

0.072 

P10 

10.4945 

0.410 

0.129 

0.386 

0.145 

0.393 

0.140 

0.609 

0.048 

P14 

10.5321 

0.417 

0.124 

0.400 

0.135 

0.404 

0.133 

0.620 

0.045 

P18 

10.5711 

0.432 

0.115 

0.  392 

0.141 

0.428 

0.118 

0.648 

0.039 

P20 

10.5911 

0.440 

0.111 

0.400 

0.135 

0.398 

0.137 

0.622 

0.045 

P24 

10.6321 

0.430 

0.116 

0.385 

0.146 

0.395 

0.139 

0.616 

0.046 

P28 

10.6746 

0.640 

0.041 

0.355 

0.169 

0.380 

0.150 

0.601 

0.050 

P32 

10.7186 

0.402 

0.134 

0.340 

0.183 

0.350 

0.174 

0.570 

0.058 

P36 

10.7641 

0.625 

0.044 

0.538 

0.068 

0.464 

0.098 

0.787 

0.020 

P38 

10.7874 

0.420 

0.122 

0.320 

0.202 

0.316 

0.206 

0.534 

0.069 

P40 

10.8112  s 

>1.66 

<0.0003 

0.308 

0.214 

0.527 

0.072 

P42 

10.8352 

0.455 

0.103  f 

0.350 

0.174 

P46 

10.8847 

0.463 

0.099 

0.357 

0.168 

*:  Values  scaled  from  NRL  graph. 

**:  Values  from  AFCRL-TR-72-0611 ,  with  y  -  0.013  added  for  normalized  "clear" 
atmosphere  (from  "average  continenta?  aerosol  model"). 


Section  A- 9 


TRANSMITTANCE  OF  0.6525  MICROMETER  LASER  BEAM  OVER 
66-KILOMETER  SLANT  PATH 

In  Section  A-5  the  following  equivalent  sea  level  path  lengths 
(ESLP's)  are  determined  for  atmospheric  aerosol  partial  transmittances : 

1.  Normalized  "clear"  atmosphere  (sea  level  visibility  =  23  km): 

ESLP,  0-1  km  vertical,  =  0.664  y 

3 

ESLP,  1-12  km  vertical,  =  0.656  y 

a 

2.  Normalized  "hazy"  atmosphere  (sea  level  visibility  =  5  km): 

ESLP,  0-1  km  vertical,  =  0.604  y  . 

a(h) 

ESLP,  1-12  km  vertical,  =  0.373  y 

a  (h) 

Also,  y^  (hazy)  =  4.87  y^  (clear). 

The  following  values  of  y  at  0.6525  micrometer  were  determined  by 

3 

linear  interpolation  between  tabled  values  for  0.63  micrometer  and  0.69 
micrometer,  from  data  supplied  by  AFCRL  for  their  "average  continental 
aerosol  model"  (1974)  and  "estimated  average  marine  aerosol  model"  (1975) 
in  private  communications  of  3  March  1975  (Ref.  19)  and  25  March  1975  (Ref. 20), 
respectively : 

y  (continental,  "clear")  =  0.1415  per  km 

3. 

y  (marine,  "clear")  =  0.1471  per  km 

3 

The  "optical  depth"  (OD)  of  scattering  due  to  atmospheric  gases  from 
0  to  12  km  altitude  (vertical)  was  determined  from  Ref.  1,  at  0.6328  micro¬ 
meter  and  0.6943  micrometer,  for  the  tropical,  mid-latitude  summer,  and  mid- 

latitude  winter  atmospheric  models.  Interpolation  to  0.6525  micrometer  was 

-4 

accomplished  using  the  X  dependence  of  molecular  scattering,  with  the 
following  results; 

OD,  0-12  km,  tropical:  0.0390 

OD,  0-12  km,  mid-latitude  summer:  0.0392 

OD,  0-12  km,  mid-latitude  winter  0.0402 

The  foregoing  values  and  relationships  were  used  to  develop  Tables 

A-20  and  A-21,  and  the  results  are  graphed  in  Figure  A-32. 


!sum  =  0. 977  y  ,,  , 
a  (h) 
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Table  A-20.  Continental  Aerosols 


Latitude 

/Season 

Model 

Aerosol 

Condition 

Eho 

m 

0-12  km 

Thy 

a 

0-12  km 
(=  by)* 

a 

£hy  t 
0-12  km 

OD  = 
5.41Ehyt 

T= 

-OD 

e 

T 

Clear 

0.0390 

0.1868 

0.2258 

1.2216 

0.295 

MS 

0.0392 

0.2260 

1.223 

0.294 

MW 

0.0402 

0.2270 

1.228 

0.293 

T 

Lt.Haze 

0.3547# 

0.3937 

2.130 

0.119 

MS 

0.3939 

2.131 

0.119 

MW 

0.3949 

2.136 

0.118 

T 

Hazy 

0.6735 

0.7125 

3.855 

0.021 

MS 

0.7127 

3.856 

0.021 

MW 

0.7137 

3.861 

0.021 

*:  b  =  1.32  (clear),  4.76  (hazy). 

//:  Light  haze  Ehya  =  geometric  mean  of  clear  and  hazy  values. 


Table  A-21.  Composite  Marine  Aerosols 


Latitude 

/Season 

Model 

Aerosol 

Condition 

Eha 

m 

0-12  km 

Thy  ** 
a 

0-1  km 
(marine) 

Ehy  ** 
a 

1-12  km 
(continental) 

Ehyt 

0-12  km 

OD  = 
5.41Ehyt 

T= 

-OD 

e 

T 

Clear 

0.0390 

0.0977 

0.0928 

0.2295 

1.2416 

0.289 

MS 

0.0392 

0.2297 

1.2427 

0.289 

MW 

0.0402 

0.2307 

1.2481 

0.287 

T 

Lt.Haze 

0.2056 

0.1544 

0.3990 

2.1586 

0.115 

MS 

0.3992 

2.1597 

0.115 

MW 

0.4002 

2.1651 

0.115 

T 

Hazy 

0.4327 

0.2570 

0.7287 

3.9423 

0.019 

MS 

0.7289 

3.9433 

0.019 

MW 

0.7299 

3.9488 

0.019 

**:  See 

text  for  conversion 

factors  and 

sea  level  coefficients. 
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Clear:  sea  level  visibility  =  23  km 
Light  haze:  sea  level  visibility  =  10.8  km 
Hazy:  sea  level  visibility  =  5  km 
W  =  mid-latitude  winter  model 
S  =  mid-latitude  summer  model 
T  =  tropical  model 


Average  Continental 
Atmosphere 


Estimated  Composite 
Marine  Atmosphere 


Clear 
W  S  T 


Clear 
W  S  T 


0.3 


Figure  A-32.  Calculated  Transmittance  of  0.6525  ym 
Laser  Over  66-Kilometer  Slant  Path 


Section  A-10 


CALCULATED  66-KILOMETER  SLANT  PATH  TRANSMITTANCES 
FOR  36  DF  LASER  LINES 


INTRODUCTION 

The  atmospheric  transmittance  of  a  66-kilometer  slant  path  was 
calculated  for  36  DF  laser  lines  covering  the  frequency  range  from  2419.070 
waves  per  centimeter  to  2862.653  waves  per  centimeter.  The  corresponding 
wavelength  range  is  from  3.493263  micrometers  to  4.133820  micrometers.  The 
specified  slant  path  of  interest  is  from  sea  level  to  40,000  feet  altitude, 
with  a  horizontal  range  of  35  nautical  miles.  Corresponding  metric  distances 
are  0  to  12.192  kilometers  altitude,  64.82  kilometers  horizontal  range,  and 
65.96  kilometers  slant  range,  with  path  elevation  angle  =  10.65°.  (All 
calculations  were  accomplished  in  the  metric  system.) 

RESULTS 

The  results  are  shown  in  Tables  A-22,  23,  and  24,  for  18  season, 
latitude,  and  atmospheric  aerosol  combinations. 

The  headings  of  "MW,"  "MS,"  and  "tropical"  represent  normalized 
mid-latitude  winter,  mid-latitude  summer,  and  tropical  atmospheres,  re¬ 
spectively,  as  specified  by  McClatchey,  et  al,  in  Optical  Properties  of  the 
Atmosphere  (Third  Edition),  Report  No.  AFCRL^72-0497  (Ref.  1). 

The  terms  "clear,"  "light  haze,"  and  "hazy"  in  the  table  titles 
indicate  the  levels  of  atmospheric  aerosols  added  to  the  clean  air  models, 
and  represent  normalized  aerosol  distributions  and  concentrations  for  clear  , 
moderately  hazy  and  hazy  atmospheric  conditions,  respectively.  The  cor¬ 
responding  eye  visibility  ranges  at  sea  level  are  23  km,  10.8  km,  and  5  km, 
respectively. 

The  headings  of  "continental  aerosols"  and  "composite  marine  aero¬ 
sols"  indicate  the  source  of  the  aerosol  constituents.  The  composite  marine 
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aerosol  model  incorporates  a  high  concentration  of  sea  spray  in  the  first 
kilometer  above  sea  level. 


CALCULATION  METHOD 

1.  Revised  laser  line  frequencies  were  employed,  per  advance  un¬ 
published  data  (Ref.  16)  received  from  Dr.  Ron  Long,  OSU  Electro-Science 
Laboratory  on  8  August  1975.  The  stated  accuracy  of  these  revised  frequencies 
is  +  0.003  waves  per  centimeter. 

2.  Revised  tables  of  molecular  absorption  coefficients  by  altitude 
layers  (for  the  OSU  line  frequencies)  were  used,  per  advance  unpublished 
data  (Ref.  41)  received  from  AFCRL  on  11  September  1975.  Computer  calcula¬ 
tion  of  these  tables  was  accomplished  by  AFCRL  upon  request  of  SRL  via  our 
letter  of  11  August  1975  to  Dr.  John  Selby,  Optical  Physics  Laboratory 

of  AFCRL. 


3.  Revised  tables  of  aerosol  absorption  and  scattering  coefficients 
for  an  "average  continental  aerosom  model"  were  employed,  per  unpublished 
data  (Ref.  19)  received  from  Dr.  Selby  by  letter  of  3  March  1975.  SRL  de¬ 
rived  vertical  "optical  depth"  values  (for  our  nominal  slant  path)  from 
these  revised  sea  level  tables  per  techniques  deduced  from  AFCRL-72-0497 
(Ref.  1). 


4.  Tables  of  estimated  aerosol  absorption  and  scattering  coefficients 
for  an  "estimated  marine  aerosol  model"  were  employed,  per  advance  unpublished 
data  (Ref.  20)  received  from  Dr.  Selby  by  letter  of  25  March  1975.  Again, 

SRL  derived  vertical  optical  depth  values  from  these  tables.  The  "estimated 
marine  aerosol  model"  is  composed  of  60%  sea  spray  and  40%  continental 
aerosols . 


5.  SRL's  "composite"  marine  atmosphere  assumes  that  the  AFCRL 
"estimated  marine  aerosol  model"  (60%  sea  spray)  exists  from  sea  level 
to  one  kilometer  altitude,  and  continental  aerosols  exist  above  one 
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kilometer.  Both  are  appropriately  tapered  with  altitude,  using  factors 
derived  from  AFCRL-72-0497  (Ref.  1). 

6.  Values  for  H^O  continuum  absorption  (omitted  from  the  AFCRL 
tables)  were  taken  from  revised  OSU  advance  unpublished  data  (Ref.  15) 
attributed  to  Burch.  These  sea  level  values  were  converted  to  optical 
depth  values  for  vertical  paths  by  SRL,  using  the  methods  and  data  of 
AFCRL-72-0497  (Ref.  1). 

7.  The  vertical  optical  depth  (OD)  values  (0-12  km),  computed 
for  each  class  of  contributor,  were  summed  to  determine  a  total  vertical 
OD.  The  path  OD  was  then  taken  as  5.41  times  the  total  vertical  OD,  since 
sec  (90°  -  10.65°)  is  5.41. 

NOTE: 

The  AFCRL  marine  atmosphere  aerosol  model  and  the  assumptions  of 
SRL's  composite  aerosol  model  are  estimates  based  upon  the  limitations 
of  currently  available  data.  An  extensive  at-sea  measurements  program 
would  be  needed  to  improve  these  estimates  significantly. 
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Table  A-22.  Calculated  66-Kilometer  Slant  Path 
Transmittances  for  36  DF  Laser 
Lines,  Clear  Atmosphere 


Line  ID 

Frequency 

(cm"l) 

Continental  Aerosols 

Composite  Marine 

Aerosols 

MW 

MS 

Tropical 

MW 

MS 

Tropical 

3-2  P (13) 

2419.070 

0.173 

0.132 

0.109 

0.143 

0.109 

0.090 

3-2  P ( 1 2 ) 

2445.356 

0.223 

0.172 

0.145 

0.184 

0.142 

0,120 

3-2  P(ll) 

2471.245 

0.321 

0.257 

0.218 

0.265 

0.212 

0.180 

3-2  P (10) 

2496.721 

0.507 

0.429 

0.377 

0.417 

0.353 

0.310 

2-1  P ( 1 3) 

2500.428 

0.531 

0.453 

0.400 

0.437 

0.373 

0.329 

3-2  P(9) 

2521.769 

0.638 

0.556 

0.497 

0.525 

0.457 

0.409 

2-1  P(12) 

2527.391 

0.656 

0.574 

0.515 

0.539 

0.472 

0.423 

3-2  P(8) 

2546.375 

0.448 

0.404 

0.365 

0.368 

0.332 

0.300 

2-1  P(ll) 

2553.953 

0.625 

0.563 

0.510 

0.513 

0.462 

0.419 

3-2  P(7) 

2570.522 

0.184 

0.180 

0.165 

0.151 

0.148 

0.135 

2-1  P(10) 

2580.097 

0.154 

0.142 

0.130 

0.126 

0.116 

0.106 

1-0  P (13) 

2583.486 

0.511 

0.448 

0.403 

0.419 

0.367 

0.330 

3-2  P(6) 

2594.198 

0.768 

0.658 

0.584 

0.629 

0.539 

0.478 

2-1  P (9) 

2605.807 

0.772 

0.616 

0.524 

0.631 

0.504 

0.429 

1-0  P(12) 

2611.142 

0.807 

0.699 

0.623 

0.660 

0.571 

0.509 

3-2  P(5) 

2617.386 

0.816 

0.726 

0.657 

0.667 

0.593 

0.537 

2-1  P(8) 

2631.068 

0.752 

0.601 

0.509 

0.614 

0.491 

0.415 

1-0  P(ll) 

2638.392 

0.403 

0.067 

0.023 

0.329 

0.055 

0.019 

3-2  P(4) 

2640.074 

0.761 

0.549 

0.438 

0.621 

0.448 

0.358 

2-1  P(7) 

2655.863 

0.669 

0.371 

0.255 

0.546 

0.302 

0.208 

3-2  P(3) 

2662.246 

0.793 

0.625 

0.525 

0.647 

0.509 

0.428 

1-0  P(10) 

2665.219 

0.760 

0.614 

0.523 

0.619 

0.500 

0.427 

2-1  P(6) 

2680.179 

0.747 

0.521 

0.407 

0.608 

0.424 

0.332 

3-2  P(2) 

2683.890 

0.821 

0.709 

0.625 

0.668 

0.577 

0.509 

1-0  P(9) 

2691.607 

2-1  P(5) 

2703.999 

0.827 

0.698 

0.605 

0.612 

0.567 

0.492 

1-0  P (8) 

2717.539 

0.675 

0.328 

0.205 

0.548 

0.266 

0.167 

2-1  P(4) 

2727.309 

0.771 

0.541 

0.419 

0.625 

0.439 

0.340 

1-0  P (7) 

2742.998 

0.741 

0.498 

0.309 

0.601 

0.403 

0.250 

2-1  P(3) 

2750.094 

0.781 

0.617  . 

0.511 

0.633 

0.500 

0.414 

1-0  P (6) 

2767.968 

0.517 

0.323 

0.231 

0.419 

0.261 

0.187  ■ 

2-1  P  (2) 

2772.340 

0.0563 

0.00016 

<.00001 

0.046 

0.00013 

<.00001 

1-0  P ( 5) 

2792.434 

0.728 

0.479' 

0.352 

0.588 

0.387 

0.284 

1-0  P(4) 

2816.380 

0.607 

0.311 

0.194 

0.489 

0.250 

0.156 

1-0  P(3) 

2839.791 

0.742 

0.445 

0.305 

0.596 

0.358 

0.245 

1-0  P(2) 

2862.653 

0.715 

0.407 

0.268 

0.574 

0.326 

0.215 

Note:  Calculated  values  are  omitted  for  the  1-0  P(9)  line  because  an  erroneous 
frequency  value  was  employed  by  AFCRL  in  the  computer  calculation  of  coeffi¬ 
cients  for  this  line,  making  the  results  meaningless.  Dr.  Selby  reported  (Ref. 43) 
(Ref.  43)  that  he  is  checking  H20  line  strengths  throughout  this  spectral  region 
and  will  supply  a  complete  set  of  new  tables  eventually,  in  lieu  of  re-running 
1-0  P(9) . 
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Table  A-23.  Calculated  66-Kilometer  Slant  Path 
Transmittances  for  36  DF  Laser  Lines, 
Light  Haze  Atmosphere 


Frequency  Continental  Aerosols  Composite  Marine  Aerosols 


Line  ID 

(cm-l) 

MW 

MS 

Tropical 

MW 

MS 

Tropical 

3-2  P(13) 

2419.070 

0.154 

0.118 

0.097 

0.104 

0.079 

0.065 

3-2  P(12) 

2445.356 

0.199 

0.153 

0.130 

0.133 

0.102 

0.087 

3-2  P(ll) 

2471.245 

0.286 

0.229 

0.194 

0.191 

0.153 

0.130 

3-2  P(10) 

2496.721 

0.451 

0.382 

0.335 

0.300 

0.254 

0.223 

2-1  P ( 1 3) 

2500.428 

0.472 

0.403 

0.356 

0.314 

0.268 

0.236 

3-2  P ( 9) 

2521.769 

0.567 

0.494 

0.442 

0.376 

0.328 

0.293 

2-1  P(12) 

2527.391 

0.583 

0.510 

0.457 

0.386 

0.338 

0.303 

3-2  P(8) 

2546.375 

0.398 

0.358 

0.324 

0.263 

0.237 

0.214 

2-1  P(ll) 

2553.953 

0.555 

0.500 

0.453 

0.366 

0.330 

0.299 

3-2  P(7) 

2570.522 

0.163 

0.160 

0.146 

0.108 

0.105 

0.096 

2-1  P(10) 

2580.097 

0.137 

0.126 

0.115 

0.090 

0.083 

0.076 

1-0  P(13) 

2583.486 

0.453 

0.397 

0.357 

0.298 

0.261 

0.235 

3-2  P(6) 

2594.198 

0.680 

0.583 

0.517 

0.447 

0.383 

0.340 

2-1  P ( 9 ) 

2605.807 

0.683 

0.545 

0.464 

0.448 

0.358 

0.304 

1-0  P(12) 

2611.142 

0.714 

0.618 

0.551 

0.468 

0.405 

0.361 

3-2  P(5) 

2617.386 

0.722 

0.642 

0.581 

0.473 

0.420 

0.380 

2-1  P(8) 

2631.068 

0.665 

0.532 

0.450 

0.434 

0.347 

0.294 

1-0  P(ll) 

2638.392 

0.357 

0.060 

0.020 

0.232 

0.039 

0.013 

3-2  P(4) 

2640.074 

0.672 

0.485 

0.387 

0.439 

0.316 

0.253 

2-1  P(7) 

2655.863 

0.591 

0.327 

0.225 

0.385 

0.213 

0.147 

3-2  P(3) 

2662.246 

0.700 

0.551 

0.463 

0.456 

0.359 

0.302 

1-0  P(10) 

2665.219 

0.670 

0.542 

0.462 

0.436 

0.352 

0.300 

2-1  P(6) 

2680.179 

0.658 

0.459 

0.359 

0.427 

0.298 

0.233 

3-2  P(2) 

2683.890 

0.723 

0.625 

0.551 

0.469 

0.405 

0.358 

1-0  P(9) 

2691.607 

2-3  P(5) 

2703.999 

0.728 

0.615 

0.533 

0.471 

0.398 

0.345 

1-0  P(8) 

2717.539 

0.594 

0.289 

0.181 

0.383 

0.186 

0.117 

2-1  P(4) 

2727.309 

0.680 

0.476 

0.369 

0.437 

0.307 

0.238 

1-0  P ( 7) 

2742.998 

0.652 

0.437 

0.272 

0.419 

0.281 

0.175 

2-1  P  ( 3) 

2750.094 

0.686 

0.542 

0.449 

0.441 

0.348 

0.289 

1-0  P(6) 

2767.968 

0.454 

0.284 

0.203 

0.291 

0.182 

0.130 

2-1  P  ( 2) 

2772.340 

0.049 

0.00014 

<.00001 

0.032 

0.00009 

<.00001 

1-0  P ( 5) 

2792.434 

0.639 

0.420 

0.309 

0.407 

0.268 

0.197 

1-0  P(4) 

2816.380 

0.512 

0.272 

0.170 

0.338 

0.173 

0.108 

1-0  P ( 3) 

2839.791 

0.648 

0.389 

0.267 

0.410 

0.246 

0.169 

1-0  P (2) 

2862.653 

0.625 

0.356 

0.234 

0.393 

0.223 

0.147 

See  Note  for  Table  A-22. 
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Table  A- 24.  Calculated  66-Kilometer  Slant  Path  Transmittances 
for  36  DF  Laser  Lines,  Hazy  Atmosphere 


Line  ID 


3-2  P(13) 
3-2  P(12) 
3-2  P(ll) 
3-2  P(10) 

2- 1  P(13) 

3- 2  P(9) 

2- 1  P (1 2 ) 

3- 2  P ( 8) 

2- 1  P(ll) 

3- 2  P(7) 

2- 1  P(10) 

1- 0  P (13) 

3- 2  P (6) 

2- 1  P(9) 

1- 0  P(12) 

3- 2  P(5) 

2- 1  P (8) 

1- 0  P(ll) 

3- 2  P(4) 

2- 1  P ( 7) 

3- 2  P ( 3) 

1- 0  P(10) 

2- 1  P(6) 

3- 2  P(2) 

1- 0  P (9 ) 

2- 1  P(5) 

1- 0  P (8) 

2- 1  P(4) 

1- 0  P(7) 

2- 1  P ( 3) 

1- 0  P (6 ) 

2-  1  P  (2) 
1-0  P(5) 
1-0  P(4) 
1-0  P(3) 
1-0  P(2) 


See  Note 


Frequency  Continental  Aerosols _  Composite  Marine  Aerosols 


(cm-1) 

MW 

MS 

Tropical 

MW 

MS 

Tropical 

2419.070 

0.125 

0.095 

0.078 

0.054 

0.041 

0.034 

2445.356 

0.160 

0.124 

0.105 

0.069 

0.053 

0.045 

2471.245 

0.230 

0.184 

0.156 

0.098 

0.078 

0.067 

2496.721 

0.362 

0.306 

0.269 

0.153 

0.130 

0.114 

2500.428 

0.379 

0.323 

0.285 

0.160 

0.137 

0.121 

2521.769 

0.454 

0.396 

0.353 

0.191 

0.166 

0.149 

2527.391 

0.466 

0.408 

0.366 

0.196 

0.171 

0.154 

2546.375 

0.317 

0.286 

0.259 

0.133 

0.119 

0.108 

2553.953 

0.442 

0.398 

0.361 

0.184 

0.166 

0.150 

2570.522 

0.130 

0.127 

0.116 

0.054 

0.053 

0.048 

2580.09/ 

0.109 

0.100 

0.091 

0.045 

0.042 

0.038 

2583.486 

0.360 

0.315 

0.284 

0.149 

0.130 

0.118 

2594.198 

0.541 

0.463 

0.411 

0.223 

0.191 

0.169 

2605.807 

0.542 

0.433 

0.368 

0.223 

0.178 

0.151 

2611.142 

0.566 

0.490 

0.457 

0.232 

0.201 

0.179 

2617.386 

0.572 

0.509 

0.460 

0.234 

0.208 

0.189 

2631.068 

0.527 

0.421 

0.356 

0.215 

0.172 

0.145 

2638.392 

0.282 

0.047 

0.016 

0.115 

0.019 

0.007 

2640.074 

0.532 

0.383 

0.306 

0.217 

0.156 

0.125 

2655.863 

0.467 

0.259 

0.178 

0.189 

0.105 

0.072 

2662.246 

0.553 

0.435 

0.366 

0.233 

0.178 

0.148 

2665.219 

0.329 

0.427 

0.365 

0.214 

0.173 

0.147 

2680.179 

0.519 

0.362 

0.283 

0.209 

0.146 

0.115 

2683.890 

0.570 

0.492 

0.434 

0.229 

0.198 

0.174 

2691.607 

2703.999 

0.572 

0.483 

0.419 

0.228 

0.193 

0.167 

2717.539 

0.466 

0.227 

0.142 

0.185 

0.090 

0.056 

2727.309 

0.532 

0.373 

0.289 

0.211 

0.148 

0.115 

2742.998 

0.511 

0.344 

0.213 

0.201 

0.135 

0.084 

2750.094 

0.537 

0.425 

0.352 

0.211 

0.167 

0.138 

2767.968 

0.355 

0.222 

0.158 

0.139 

0.087 

0.062 

2772.340 

0.039 

0.00011 

<.00001 

0.015 

<.00001 

<.00001 

2792.434 

0.498 

0.327 

0.241 

0.193 

0.127 

0.093 

2816.380 

0.414 

0.212 

0.132 

0.159 

0.081 

0.051 

2839.791 

0.503 

0.302 

0.207 

0.'191 

0.115 

0.079 

2862.653 

0.484 

0.275 

0.181 

0.182 

0.103 

0.068 

for  Table  A-22. 
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Section  A- 11 


CALCULATED  AVERAGE  TRANSMITTANCE  OF  A  500-WATT  DF  LASER 
OVER  66-KILOMETER  SLANT  PATH 

The  attached  Tables  A-25,  26,  27  and  28,  and  Figure  A-33  show  the 
results  of  calculations  of  the  transmittance  of  the  output  beam  of  the  500- 
watt  HAC/AFAL  laser.  The  power  distribution  of  that  laser  was  obtained 
from  a  data  sheet  (Ref.  44)  supplied  by  Dr.  Ron  Paulson,  AFAL. 

Revised  laser  line  frequencies  were  employed,  per  advance  unpublished 
data  (Ref.  16)  received  from  Dr.  Ron  Long,  OSU  Electro-Science  Laboratory  on 
8  August  1975. 

Revised  tables  of  molecular  absorption  coefficients  versus  altitude 
(for  the  OSU  line  frequencies)  were  used,  per  advance  unpublished  data 
(Ref.  41)  received  from  AFCRL  on  11  September  1975. 

Revised  tables  of  aerosol  absorption  and  scattering  coefficients 
for  an  "average  continental  aerosol  model"  were  employed,  per  unpublished 
data  (Ref.  19)  received  from  Dr.  Selby  by  letter  of  3  March  1975.  SRL 
derived  vertical  "optical  depth"  values  (for  our  nominal  slant  path)  from 
these  revised  tables  per  techniques  deduced  from  AFCRL-72-0497  (Ref.  1). 

Tables  of  estimated  aerosol  absorption  and  scattering  coefficients 
for  an  "estimated  marine  aerosol  model"  were  employed,  per  advance  unpub¬ 
lished  data  (Ref.  20)  received  from  Dr.  Selby  by  letter  of  25  March  1975. 
Again,  vertical  optical  depth  values  were  derived  from  these  tables.  The 
"estimated  marine  aerosol  model"  is  composed  of  60%  sea  spray  and  40% 
continental  aerosols. 

SRL's  "composite"  marine  atmosphere  assumes  that  the  AFCRL  "estimated 
marine  aerosol  model"  (60%  sea  spray)  exists  from  sea  level  to  one  kilometer 
altitude,  and  continental  aerosols  exist  above  one  kilometer.  Both  are 
appropriately  tapered  with  altitude,  using  factors  derived  from  AFCRL-72- 
0497  (Ref.  1). 
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"Clear"  and  "hazy"  on  the  tables  and  graph  refer  to  a  normalized 
clear  atmosphere  (visibility  =  23  kilometers)  and  normalized  hazy  atmos¬ 
phere  (visibility  =  5  kilometers),  respectively,  at  sea  level.  "W"  and 
"MW"  represent  McClatchey's  mid-latitude  winter  atmosphere  model,  "S"  and 
"MS"  represent  the  mid-latitude  summer  model,  and  "T"  represents  the  tropical 
model,  all  per  AFCRL-72-0497  (Ref.  1). 

The  "average"  or  "total"  transmittance  values  calculated  are 
weighted  averages:  the  calculated  total  transmittance  at  each  laser  line 
frequency  is  weighted  by  the  transmitted  power  percentage  at  that  line. 

Hence  the  resulting  values  represent  the  fraction  of  the  total  transmitted 
power  surviving  at  the  end  of  the  path. 

A  word  of  caution:  The  AFCRL  marine  atmosphere  aerosol  model  and 
the  assumptions  of  SRL's  composite  aerosol  model  are  estimates  based  upon 
the  limitations  of  currently  available  data.  An  extensive  at-sea  measure¬ 
ments  program  would  be  needed  to  improve  these  estimates  significantly. 
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Table  A-25.  Average  Transmittance  Calculation  for  500-Watt 
DF  Laser  over  66-Kilometer  Slant  Path: 

Clear  Atmosphere,  Continental  Aerosols 


It 


Model  Atmosphere 


Line 

%* 

MW  - 

Clear 

MS  - 

Clear 

T  - 

Clear 

Id 

PT 

T 

/  p  * 

R 

T 

%PR* 

T 

/o  rR 

1-0  P (4) 

0.3 

0.607 

0.18 

0.311 

0.09 

0.194 

0.06 

1-0  P (5) 

8.38 

0.728 

6.10 

0.479 

4.01 

0.352 

2.95 

1-0  P (6) 

9.68 

0.517 

5.00 

0.323 

3.13 

0.231 

2.24 

1-0  P (7) 

8.58 

0.741 

6.36 

0.498 

4.27 

0.  309 

2.65 

2-1  P (4) 

0.70 

0.771 

0.54 

0.541 

0.38 

0.419 

0.29 

1-0  P (8) 

3.99 

0.675 

2.69 

0.328 

1.31 

0.205 

0.82 

2-1  P (5) 

8.68 

0.827 

7.18 

0.698 

6.06 

0.605 

5.25 

2-1  P (6) 

9.98 

0.747 

7.46 

0.521 

5.20 

0.407 

4.06 

2-1  P (7) 

9.48 

0.669 

6.34 

0.371 

3.52 

0.255 

2.42 

2-1  P(8) 

8.68 

0.752 

6.53 

0.601 

5.22 

0.509 

4.42 

3-2  P (5) 

8.08 

0.816 

6.  ife 

0.726 

5.87 

0.657 

5.31 

2-1  P (9) 

1.0 

0.772 

0.77 

0.616 

0.62 

0.524 

0.52 

3-2  P (6) 

9.08 

0.768 

6.97 

0.658 

5.97 

0.584 

5.30 

3-2  P(7) 

8.48 

0.184 

1.56 

0.180 

1.53 

0.165 

1.40 

3-2  P (8) 

4.69 

0.448 

2.10 

0.404 

1.89 

0.365 

1.71 

3-2  P (9) 

0.2 

0.638 

0.13 

0.556 

0.11 

0.497 

0.10 

Totals : 

100.0% 

66.5% 

49.2% 

39.5% 

*:  P_  and 
T 

P  expressed  as  % 
K 

of  total 

power  transmitted. 

P,r  is 

trans- 

mitted  power;  P  is  beam  power  at  end  of  path. 
K 
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Table  A-26.  Average  Transmittance  Calculation  for  500-Watt 
DF  Laser  over  66-Kilometer  Slant  Path: 

Hazy  Atmosphere,  Continental  Aerosols 


Line 

%* 

Model  Atmosphere 

MW  - 

Hazy 

MS  - 

Hazy 

T  - 

Hazy 

Id 

PT 

T 

% 

_ 

V 

T 

%  V 

T 

/  p  * 
R 

1-0 

P(4) 

0.3 

0.414 

0. 

12 

0.212 

0.06 

0.132 

0.04 

1-0 

P(5) 

8.38 

0.498 

4. 

17 

0.327 

2.74 

0.241 

2.02 

1-0 

P  (6) 

9.68 

0.355 

3. 

44 

0.222 

2.15 

0.158 

1.53 

1-0 

P(7) 

8.58 

0.511 

4. 

38 

0.  344 

2.95 

0.213 

1.83 

2-1 

P(4) 

0.  70 

0.532 

0. 

37 

0.  373 

0.26 

0.289 

0.20 

1-0 

P(8) 

3.99 

0.466 

1. 

86 

0.227 

0  .91 

0.142 

0.57 

2-1 

P  (5) 

8.68 

0.572 

4. 

96 

0.483 

4.19 

0.419 

3.64 

2-1 

P  (6) 

9.98 

0.519 

5. 

18 

0.362 

3.61 

0.283 

2.82 

2-1 

P(7) 

9.48 

0.467 

4. 

43 

0.259 

2.46 

0.178 

1.69 

2-1 

P(8) 

8.68 

0.527 

4. 

57 

0.421 

3.65 

0.356 

3.09 

3-2 

P(5) 

8.08 

0.572 

4 . 

62 

0.509 

4.11 

0.460 

3.72 

2-1 

P(9) 

1.0 

0.542 

0. 

54 

0.433 

0.43 

0.368 

0.37 

3-2 

P(6) 

9.08 

0.541 

4. 

91 

0.463 

4.20 

0.411 

3.73 

CN 

1 

P(7) 

8.48 

0.130 

1. 

10 

0.127 

1.08 

0.116 

0.98 

3-2 

P  (8) 

4.69 

0.317 

1. 

49 

0.286 

1.34 

0.259 

1.21 

3-2 

P(9) 

0.2 

0.454 

0. 

09 

0.396 

0.08 

0.353 

0.07 

Totals: 

100. 0% 

46. 

3% 

34.2% 

27.5% 

*:  1 

and 

T 

P  expressed  as  % 

of  total 

power  transmitted. 

P  is  trans- 

mitted  power;  P  is  beam  power  at  end  of  path. 
K 


Table  A-27.  Average  Transmittance  Calculation  for  500-Watt 
DF  Laser  over  66-Kilometer  Slant  Path: 

Clear  Atmosphere,  Composite  Aerosols 


Model  Atmosphere 


L:fne 

Id 

%* 

P 

MW  - 

Clear 

MS  - 

Clear 

T  - 

Clear 

T 

%PR* 

T 

1  %  P  * 

1  *R 

T 

1-0  P (4) 

0.3 

0.489 

0.15 

0.250 

0.08 

0.156 

0.05 

1-0  P (5) 

8.38 

0.  588 

4.93 

0.387 

3.24 

0.284 

2.38 

1-0  P (6) 

9.68 

0.419 

4.06 

0.261 

2.53 

0.187 

1.81 

1-0  P(7) 

8.58 

0.601 

5.16 

0.403 

3.46 

0.250 

2.15 

2-1  P  (4) 

0.70 

0.625 

0.44 

0.439 

0.31 

0.340 

0.24 

1-0  P (8) 

3.99 

0.548 

2.19 

0.266 

1.06 

0.167 

0.67 

2-1  P (5) 

8.68 

0.612 

5.31 

0.567 

4.92 

0.492 

4.27 

2-1  P (6) 

9.98 

0.608 

6.07 

0.424 

4.23 

0.332 

3.31 

2-1  P ( 7 ) 

9.48 

0.546 

5.18 

0.302 

2.86 

0.208 

1.97 

2-1  P (8) 

8.68 

0.614 

5.33 

0.491 

4.26 

0.415 

3.60 

3-2  P (5) 

8.08 

0.667 

5.39 

0.593 

4.79 

0.537 

4.34 

2-1  P (9) 

1.0 

0.631 

0.63 

0.504 

0.50 

0.429 

0.43 

3-2  P (6) 

9.08 

0.629 

5.71 

0.539 

4.89 

0.478 

4.34 

3-2  P(7) 

8.48 

0.151 

1.28 

0.148 

1.26 

0.135 

1.14 

3-2  P (8) 

4.69 

0.368 

1.73 

0.332 

1.56 

0.300 

1.40 

3-2  P (9) 

0.2 

0.525 

0.11 

0.457 

0.09 

0.409 

0.08 

Totals : 

100.0% 

53.6% 

40.0% 

32.1% 

*:  P  and  P  expressed  as  %  of  total  power  transmitted.  PT  is  trans- 
T  R 

mitted  power;  P  is  beam  power  at  end  of  path. 

K 
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Table  A-28.  Average  Transmittance  Calculation  for  500-Watt 
DF  Laser  over  66-Kilometer  Slant  Path: 

Hazy  Atmosphere,  Composite  Aerosols 


Model  Atmosphere 


Line 

%* 

MW 

-  Hazy 

MS 

-  Hazy 

T  - 

Hazy 

ID 

PT 

T 

%  P  * 

T 

%  P  * 

R 

T 

%  PR* 

1-0  P (4) 

0.  3 

0.159 

0.05 

0.081 

0.02 

0.051 

0.02 

1-0  P (5) 

8.38 

0. 193 

1.62 

0.127 

1.06 

0.093 

0.78 

1-0  P (6) 

9.68 

0.139 

1.35 

0.087 

0.84 

0.062 

0.60 

1-0  P (7) 

8.58 

0.201 

1.72 

0.220 

1.89 

0.084 

0.070 

2-1  P (4) 

0.70 

0.211 

0.15 

0.148 

0.10 

0.115 

0.08 

1-0  P (8) 

3.99 

0.185 

0.74 

0.090 

0.36 

0.056 

0.22 

2-1  P(5) 

8.68 

0.228 

1.98 

0.193 

1.68 

0.167 

1.45 

2-1  P (6) 

9.98 

0.209 

2.09 

0.146 

1.46 

0.115 

1.15 

2-1  P  (7) 

9.48 

0.189 

1.79 

0.105 

1.00 

0.072 

0.68 

2-1  P(8) 

8.68 

0.215 

1.87 

0.172 

1.49 

0.145 

1.26 

3-2  P (5) 

8.08 

0.234 

1.89 

0.208 

1.68 

0.189 

1.53 

2-1  P (9) 

1.0 

0.223 

0.22 

0.178 

0.18 

0.151 

0.15 

3-2  P (6) 

9.08 

0.223 

2.02 

0.191 

1.73 

0.169 

1.53 

3-2  P (7) 

8.48 

0.54 

0.46 

0.053 

0.45 

0.048 

0.41 

3-2  P(8) 

4.69 

0.133 

0.62 

0.119 

0.56 

0.108 

0.51 

3-2  P (9) 

0.2 

0.191 

0.04 

0.166 

0.03 

0.149 

0.03 

Totals : 

100.0% 

18.6% 

14.5% 

11.1% 

*:  PT  and 

P  expressed  as  7 
K 

of  total 

power  transmitted 

P^  is  trans- 

mitted  power;  PD  is  beam  power  at  end  of  path. 
K 
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Section  A-12 


DF  LASER  LINE-SCANNER  ATMOSPHERIC 
BACKSCATTER  PRELIMINARY  ANALYSIS 


INTRODUCTION 

A  rigorous  analytical  treatment  of  the  atmospheric  bac.kscatter  problem 
for  an  active  system  is  difficult.  For  a  vertical  or  slant  path  the  problem 
is  compounded  by  variation  in  atmospheric  scattering  and  attenuation  coeffi¬ 
cients  with  altitude. 

Problem  simplification,  including  introduction  of  a  layer-by-layer 
evaluation  and  summation  approach  as  employed  in  Section  A-3  for  path  thermal 
radiance,  offers  a  method  for  making  engineering  estimates.  This  section 
describes  such  a  method  and  gives  the  results  of  a  preliminary  "quick  look" 
analysis  whose  purpose  was  to  determine  whether  atmospheric  backscatter  needs 
to  be  considered  in  design  analysis  for  a  DF  laser  line  scanner. 

APPROACH  (see  Figure  A-34) 

Divide  the  air-to-ground  path  into  layers  for  which  the  average 
scattering  coefficients  and  attenuation  coefficients  can  be  estimated.  The 
most  convenient  layering  is  by  one-kilometer  altitude  increments  for  which 
coefficient  tables  are  available. 

Then  for  layers  other  than  the  first  (nearest  the  transmitter- 

receiver)  the  backscattered  radiance  of  the  layer  is  given  approximately 

by  the  product  of  the  incident  power,  the  path  length  in  the  layer,  and 

the  normalized  scattering  phase  function  for  180°  (at  the  wavelength  of 

•interest),  divided  by  the  nominal  beam  area  at  the  top  of  the  layer.  This 

radiance  value  is  attenuated  by  the  intervening  atmospheric  path  between 

the  layer  and  the  receiver.  That  is,  the  backscattered  radiance  from  the 
t  H 

i  layer  is  given  approximately  by 
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Figure  A- 34.  Schematic  Representation  of  Approximate  Beam/IFOV  Overlap 


.  _  P  (arl>  °1  F(180°) 

s(4V  i  rV 

4 


-TyALr.) 


W/m  -ster. 


(A-32) 


where  P  is  transmitted  power,  Ar^  is  the  path  length  in  the  i 


th 


o  and 


y,  are  the  scattering  and  attenuation  coefficients  of  the  i 

th 


layer , 
.  th 


layer,  r  is  the  total  path  length  to  the  i  layer,  6  is  the  nominal 

beam  diameter,  and  F/10no.  is  the  backscatter  normalized  scattering 

v.ioU  ; 

phase  function.  The  exponential  factor  is  squared,  to  account  for  two- 
way  atmospheric  attenuation. 


Equation  (A-32)  becomes  increasingly  inaccurate  for  atmospheric 
layers  close  to  the  source,  and  fails  completely  at  close  ranges.  Hence 
another  approach  was  employed  to  approximate  the  atmospheric  backscatter 
from  the  top  atmospheric  layer.  The  analytical  basis  for  that  approxima¬ 
tion  follows. 


Assume,  for  the  case  of  a  CW  laser  illuminator,  that  separate 
receiving  and  transmitting  apertures  are  employed  (as  would  likely  be 
necessary) ,  separated  by  a  small  amount  so  that  the  receiver  FOV  and 
transmitter  beam  path  do  not  intersect  for  at  least  the  first  few  tens  of 
meters,  say  to  range  R^.  Assume  complete  overlap  beyond  R^,  to  the 
bottom  of  the  layer,  at  range  R£.  Then,  in  terms  of  the  illustration 
below,  the  problem  is  to  write  an  expression  for  the  scattered  power  level 
P^dr  from  an  elemental  length  of  path  dr  (expressed  in  aperture  irradiance 
terms),  then  integrate  that  expression  between  limits  R^  and  R^.  The 
differential  expression  is  (see  following  page) 


Transmitter  r  dr 
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dH 

where  PD 

D 


P,  dr 
_h _ 

2 

r 


watts/meter 


2 


Po(r)Fa8o”)(!i‘Tr) 


(of  aperture) 


backscattered  power  per  unit  length  of 
beam  path 


The  scattering  coefficient  G(r)  may  t>e  assumed  constant  for  the  layer,  and 
the  attenuation  negligible.  Then  Equation  (A-32) ,  with  the  simplified  ex¬ 
pression  for  P^  inserted,  reduces  to 


dH 


PoF(180°) 


dr  ,  2 

— 2  watts/m 

r 


Integrating,  and  inserting  the  cited  limits 


H  =  PctF(180°) 


{2 


-  —  watts/m  of  aperture, 

lx-  K. 


(A-33) 


where  a  is  in  units  of  m  and  R^  and  are  in  meters. 


The  above  formula  is  undoubtedly  pessimistic  in  that  it  assumes 
maximum  overlap  of  the  receiver  IFOV  and  the  transmit  beam  occurs  abruptly 
rather  than  gradually.  On  the  other  hand,  it  ignores  backscattered  energy 
from  port  scattering  and  the  beam  sidelobes  and/or  skirts.  Because  trans¬ 
mitted  energy  from  these  sources  can  intersect  the  receiver  IFOV  very  close 
to  the  apertures,  it  appears  possible  that  these  could  be  the  strongest 
sources  of  backscattered  irradiance  of  the  receiving  aperture. 

EXAMPLE  APPLICATION 
Assumptions 

Assume  a  wavelength  of  3.73  micrometer  (the  DF  laser  2-1  P(6) 
line)  and  a  laser  effective  beam  power  of  500  watts.  Also,  assume  a  nominal 
beam  diameter  of  66  microradians  (4.36  meters  diameter  at  target  range  of 
66  kilometers),  and  receiver  IFOV  of  30  microradians  (1.98  meters  diameter 
at  target). 
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Example  applicable  coefficients,  derived  fromAFCRL's  "average 
continental  aerosol  model"  and  AFCRh-72-0497  (Ref.  1),  are  as  follows: 


Continental 

Aerosols 


a 

a 

(0  - 

1  km,  clear) 

0.0120/km 

a 

a 

(0  - 

1  km,  hazy) 

0.0533 

0 

a 

(11 

-  12  km,  clear  or  hazy) 

0.000348 

0 

m 

(0  - 

1  km) 

1 

o 

f—i 

X 

00 

<r 

0 

m 

(11 

-  12  km) 

1.47  x  10' 

Est.  Marine 
Aerosols 

0.049/km 

0.218 

0.0014 


Aerosol  scattering:  F  (180°,  3.73  micrometer)  =  0.014 

Molecular  scattering:  F  (180°,  3.73  micrometer)  =  0.12 


It  is  clear  from  the  above  that  molecular  scattering  is  negligible 
compared  to  aerosol  scattering,  so  it  will  be  ignored  in  the  following 
examples  (for  66-km  slant  path. 

Near  Scattering  (from  top  1-km  layer) 

Assume  =  100  meters,  and  marine  aerosols,  using  Equation  (A-33): 

H  «  500  (0.0014  x  10"3K0. 014)  (y—  -  -^q) 

—8  2 

=  9.8  x  10  watts/meter  (of  aperture) 

This  is  about  two  orders  of  magnitude  higher  than  the  expected  aperture 
irradiance  from  target  reflected  energy;  hence  it  appears  that  backscatter 
from  regions  close  to  the  source/receiver  must  be  dealt  with  in  system 
design. 

Far  Scattering  [from  bottom  1-km  layer,  for  example,  with  marine  aerosols 
and  a  "hazy"  atmosphere,  using  Equation  A-32) ] : 
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500(5.41) (0.218) (0.014) (0.5) ' 
v  ..  .  .  2 
4^-4) 


=  0.136  watts/m  -ster. 


Thiv.  produces  an  aperture  irradiance  of 


0.136 


(60  x  103)2 


=  3.8  x  id"11  watts/m2 


This  value  is  somewhat  smaller  than  expected  signal  irradiance  of  the  aper¬ 
ture,  but  large  enough  to  warrant  consideration  in  detailed  system  analysis. 

Intervening  layers  will  contribute  lesser  amounts  to  the  total  back- 
scattered  irradiance  of  the  receiver,  and  should  be  considered  in  a  more 
detailed  and  comprehensive  system  analysis. 

CONCLUSION 

A  preliminary  conclusion  to  be  drawn  from  this  very  sketchy  analysis 
is  that  backscatter  from  regions  near  the  transmitter/receiver  is  apt  to  be 
most  troublesome  to  a  CW  active  scanning  system.  However,  the  signal  will 
be  largely  ac  and  the  backscatter  largely  dc,  and  this  is  one  factor  that 
can  be  used  to  discriminate  against  the  backscatter.  (Atmospheric  turbulence 
could  alter  this  discriminant  considerably.) 

Projecting  the  analysis  to  a  pulse  system,  it  appears  that  range 
gating,  to  discard  scattered  energy  from  most  of  the  path  length,  would  be 
necessary  or  desirable.  The  ac/dc  discrimination  cited  above  is,  of 
course,  not  applicable. 
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